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INTRODUCTION 


This  report  describes  investigations  of  the  interaction  of  high  power  laser 
radiationwithsurfaces  using  a  time-of-flight  mass  spectrometer  and  a 
quadrupole  mass  spectrometer.  The  work  described  in  this  report  builds 
on  and  extends  results  already  reported  in  the  semi  annual  report1'  on  this 
same  project.  Since  some  phases  of  the  investigations  are  already  completely 
reported  there,  we  will  not  duplicate  that  material.  For  background, 
however,  seme  of  the  earlier  developments  upon  which  this  report  depends 
are  summarized  below. 

Earlier  work  involved  measurements  of  the  ions  emitted  from  tungsten  targets 

and  was  carried  out  with  a  time-of-flight  mass  spectrometer.  The  species 
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emitted  at  a  laser  power  per  unit  area  of  the  order  of  50  megawatts /cm 
included  sodium,  potassium,  and  tungsten.  The  ions  were  found  to  possess 
energies  up  to  200  ev  per  ion. 

Neutral  gas  desorption  was  observed  earlier  using  a  quadrupole  mass  spectrom¬ 
eter.  The  dominant  species  found  were  adsorbed  gases  such  as  CO  and  C00, 
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in  quantities  of  the  order  of  10  molecules  per  pulse.  This  amount,  and  also 
the  manner  in  which  the  number  of  desorbed  molecules  varies  with  laser 
power,  agrees  well  with  a  thermal  interpretation  of  this  phenomenon  in  which 
the  heating  of  the  surface  by  the  laser  beam  drives  off  adsorbed  gases. 
Theoretical  work  on  possible  heating  of  the  blow-off  material,  because  of 
absorption  of  laser  light  in  inverse  Bremsstrahlung  processes  in  the  field  of 
the  ionized  atoms,  indicated  heating  is  probable  at  high  particle  densities, 
but  it  does  not  occur  at  particle  densities  representative  of  our  experimental 
conditions. 
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In  this  report  we  describe  continued  work  on  the  two  spectrometers.  We  have 
employ  d  targets  of  different  materials,  and  obtained  results  generally 
similar  to  those  found  with  tungsten  targets.  The  two  largest  topics  in  this 
report  involve  measurements  on  high  energy  neutral  particles  in  the  quadrupole 
spectrometer,  and  analysis  of  the  pulse  shapes  of  the  ions  emitted  in  the 
time-of-flight  spectrometer,  in  order  to  infer  the  velocity  distribution  of  the 
original  material.  These  two  topics  have  occupied  much  of  the  present 
report  period  and,  correspondingly,  fill  much  of  this  report. 

Section  I  discusses  the  high  speed  neutral  particle  emission.  These  pulses 
appear  rapid  compared  to  the  longer  and  slower  pulses  from  the  gas  desorption 
work.  We  shall  consider  the  various  experiments  designed  to  establish  the 
nature  of  this  emission.  We  show  that  the  surface  interaction  produces 
ultraviolet  radiation  at  wavelengths  down  to  1100A°  ,  and  high  speed  neutral 
atoms  or  molecules  with  energies  of  the  order  of  109ev.  Definite  identifications 
of  the  atomic  or  molecular  species  have  not  yet  been  attained. 

In  Section  II  we  describe  the  experiments  done  on  the  time-of-flight  spectromete 
The  main  items  here  are  investigations  on  platinum  targets,  and  an  experiment 
in  which  time  of  flight  of  the  ions  was  measured  as  a  function  of  accelerating 
voltage  as  an  additional  check  on  the  validity  of  the  theory  of  operation  of  the 
spectrometer. 

In  Section  III  we  analyze  tne  pulse  shapes  obtained  with  the  time-of-flight 
spectrometer,  and  obtain  distributions  that  indicate  a  directed  component  of 
the  ion  energy,  plus  a  smaller  random  component  superposed  on  the  directed 
component. 


Typical  results  indicate  that  the  energies  are  of  the  same  order  of  magnitude, 
but  not  equal,  for  ions  of  different  species  emitted  from  a  given  target.  The 
shapes  of  the  velocity  distributions  indicate  departures  from  a  Maxwellian 
distribution.  The  computational  procedure  used  here  rests  on  an  analysis 
developed  by  Frank  J.  Allen  of  Ballistic  Research  Laboratories  which  we 
have  programmed  for  a  computer  and  used  to  obtain  velocity  distributions 
from  the  ion  pulse.  We  include  Allen's  results  as  appendix  I  in  order  to  make 
this  report  self-contained.  Much  of  the  discussion  in  Section  III  depends  on 
Allen's  report  which  is  not  generally  available.  In  the  appendix  we  have 
condensed  Allen's  original  document,  retaining  only  those  portions  which  have 
the  most  relevance  for  this  report. 

Finally,  as  Appendix  II  we  include  a  copy  of  a  paper  based  on  work  performed 
under  this  contract  and  presented  at  the  Fourth  International  Quantum  Electronics 
Conference  in  Phoenix,  Arizona,  April  12-15,  19G6.  It  draws  mainly  on  results 
presented  previously  in  the  semi-annual  report*.  We  include  it  here  because 
it  forms  a  concise  and  unified  statement  of  work  that  has  been  pci  formed  unde*- 
this  contract. 


-3- 


I.  WORK  ON  QUADRUPOLE  MASS  SPECTROMETER 


This  section  discusses  work  done  on  the  quadrupole  mass  spectrometer 
during  the  report  period. 

The  main  work  on  this  instrument  involved  the  study  of  high  energy  molecules 
emitted  from  metallic  targets  upon  illumination  by  the  laser. 

A.  OBSERVATION  OF  THE  HIGH  SPEED  NEUTRAL  MOLECULES 

During  our  study  of  the  gas  desorption  reported  earlier,  a  pulse  was  noted 
occurring  simultaneously  with  the  laser  event.  Since  these  data  were  taken  with 
a  very  slow  time  sweep,  it  was  not  possible  to  determine  if  the  pulses  were  as 
long  as  several  tens  of  milliseconds  or  shorter.  Since  UV  light  and  high  speed 
ions  had  already  been  observed  in  the  time-of-flight  spectrometer,  it  seemed 
reasonable  that  the  pc's  observed  was  due  to  either  of  those  effects.  A 
study  was  made  to  determine  which  it  was.  When  the  pulse  was  observed  in 
higher  time  detail,  it  had  an  appearance  typified  by  the  data  shown  in  figure  1-1. 
(In  this  figure,  time  progresses  to  the  right  at  10  microseconds  per  centimeter, 
and  the  negative  going  signals  extend  upwards  with  a  sensitivity  of  10  amps 
per  cm.  )  The  sweep  is  triggered  by  the  laser  event;  not  by  the  first  pulse  of 
the  train.  The  first  pulse  is  quite  high  and  short,  beginning  with  the  laser 
event,  and  has  a  width  consistent  with  the  time  constant  of  the  detection 
system  (R  3  200  ohms  and  C  =  50  pf).  The  second  pulse  begins  to  rise 
5  usee  after  the  laser  event  and  has  a  width  of  about  7  usee.  Later  and  wider 
pulses  also  appear.  When  these  data  were  observed,  several  experiments 
were  performed  to  determine  if  the  pulses  resulted  from  ultraviolet  light, 
high  speed  ions,  or  high  speed  atoms. 


Figure  1-1  -  Photon  Pulse  and  High  Speed  Molecule  Pulse  from  a  Tungsten 
Target.  Time  progresses  to  the  right  at  10  microseconds /cm. 
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B. 


DETECTION  MECHANISM  FOR  HIGH  SPEED  MOLECULES 


Photons,  high  speed  atoms,  and  high  speed  ions  all  arrive  at  the  detector 
much  sooner  than  do  ions  resulting  from  room  temperature  gases  generated 
by  the  laser  beam.  Photons,  if  they  have  a  high  enough  energy,  and  ions  of 
nearly  any  velocity  can  be  detected  directly  at  the  multiplier.  High  speed 
neutrals,  on  the  other  hand,  pose  a  special  problem. 

When  one  is  working  with  gas  phase  neutral  molecules  moving  with  speeds 
corresponding  to  thermal  motion  at  room  temperature,  the  detection  mechanism 
usually  involves  the  ionization  of  the  atom  or  molecule  followed  by  the  sub¬ 
sequent  mass  analysis  and  detection  of  the  ion.  Under  these  conditions  the 
density  of  the  atoms  in  the  gas  phase  can  be  inferred  from  the  known  ionization 
probability  by  electron  impact  and  the  electron  beam  current.  About  1  atom 

g 

in  10  generated  is  usually  detected. 

When  high  speed  atoms  or  molecules  are  encountered,  the  same  techniques 
should  work  in  principle.  However,  the  speed  of  the  atoms,  and  the  fact  that 
they  reside  in  the  ionizing  region  a  very  short  time,  make  the  number  of  ions 
formed  by  electron  impact  very  small  compared  even  to  that  for  room 
temperature  conditions.  Because  of  the  enormous  kinetic  energy  of  the  high 
speed  atoms,  the  time  spent  in  the  ionization  region  is  very  small  (inversely 
proportional  to  velocity).  Furthermore,  the  kinetic  energy  is  dissipated 
on  the  very  first  collison  with  a  wall.  Thus  the  probability  of  detecting  a 
high  speed  atom  or  molecule  by  this  technique  is  many  orders  of  magnitude 
lower  than  that  of  detecting  thermal  atoms  and  molecules. 
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For  atoms  of  high  energy,  however,  there  is  at  least  one  detection  mechanism 
which  depends  in  no  way  on  ionization  by  an  electron  beam.  In  this  scheme 
the  atoms  are  allowed  to  strike  a  surface  and  produce  secondary  electrons. 

The  secondaries  are  then  detected  by  conventional  means.  In  our  system  the 
atoms  are  allowed  to  strike  the  first  dynode  of  the  electron  multiplier 
detector,  and  the  resulting  secondary  electrons  traverse  the  multiplier,  build¬ 
ing  up  into  an  easily  detectable  signal. 

The  number  of  atoms  in  the  initial  beam  can  be  computed  if  ''me  knows  the 

gain  of  the  electron  multiplier  and  the  value  of  Y  ,  the  coefficient  for  secondary 

electron  production  by  the  neutral  atoms  at  the  target  surface.  There  are  very 

2 

few  estimates  of  Y  in  the  literature.  One  by  Rostagni  is  reproduced  in 
table  1-1.  In  it  he  shows  a  secondary  electron  coefficient  of  0.  07  for  argon 
atoms  of  200  ev  energy.  For  neon  and  helium  he  estimates  0.  04  and  0.  39 
respectively.  Estimates  or  calculations  for  other  gases  or  for  specific 
surfaces  are  not  available.  Also  reproduced  as  part  of  table  1-1  are  data 
from  the  same  source  on  Y  for  ion  impact.  Note  how  the  values  of  Y  are 
essentially  identical  for  values  of  kinetic  energy  greater  than  200  ev.  Thus, 
for  energies  greater  than  200  ev,  it  should  be  rossible  to  infer  the  atomic 
secondary  effects  from  the  value  of  y  for  ions  on  surfaces  (At  high  kinetic 
energies  the  charge  of  the  particle  should  be  of  small  importance  when 
collision  parameters  are  involved).  Brown  gives  a  good  summary  of 
values  of  Y  in  the  literature  for  ions  on  various  surfaces.  He  indicates  that 
values  of  Y  of  the  order  of  0.  02  through  0.  2  for  energies  greater  than  200  ev 
are  the  rule. 

Since  we  will  be  discussing  the  interaction  of  hydrogen  atoms  or  molecules 
with  beryllium -copper  dynodes,  it  will  bn  necessary  .o  estimate  the  variation 
of  Y  with  atomic  or  molecular  kinetic  energy.  Seme-  sort  cf  picture  such 
as  the  following  may  b<  suggested:  The  vaiuo  of  Y  for  energy  less  than  some 
value  (perhaps  200  ev)  is  linear  with  energy.  Above  that  energy  Y  is  constant. 
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TABLE  I-  1 


VALUES  OF  THE  SECONDARY  ELECTRON  COEFFICIENT  y 


E  (ev) 


■ 

600 

400 

200 

100 

50 

30 

25 

— 

20 

Ar 

23.2 

17.0 

7.0 

1.3 

0.4 

0.02 

0.008 

Ne 

MB 

3.6 

1.0 

0.3 

0.02 

He 

m 

39.3 

8.1 

0.8 

0.06 

0.008 
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This  picture  is  consistent  with  the  assertion  that  the  kinetic  energy  of  the 
atom  or  ion  is  the  major  fraction  of  the  excitation  if  it  is  Targe  compared  to 
the  work  function  of  the  surface,  and  that  the  excitation  by  an  uncharged 
particle  is  zero  at  zero  kinetic  energy.  When  we  wish  to  interpret  the 
pulses  of  high  speed  atoms  in  relative  numbers,  we  shall  have  to  use  such 
a  picture;  however,  for  this  report  we  shall  state  only  that  this  problem  of 
the  variation  of  Ywith  atomic  kinetic  energy  is  a  serious  one. 

C.  USE  OF  THE  QUADRUPOLE  MASS  ANALYZER  TO  DETERMINE  THE 
NATURE  OF  THE  HIGH  SPEED  PULSE  TRAIN 

The  quadrupole  mass  spectrometer  is  ideally  suited  to  determine  the  nature 

of  the  train  of  pulses  observed  in  this  experiment.  It  has  been  described 
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in  a  previous  report  .  Since  it  is  a  spectrometer  with  a  line -of -sight  path 
between  the  ion  source  and  the  detector,  and  since  the  detector  is  an 
electron  multiplier  capable  of  detecting  any  particle  creating  secondary 
electrons,  it  will  allow  high  speed  ions,  excited  or  high  speed  atoms,  and 
photons  to  travel  from  the  target  to  the  detector  and  be  observed.  In  order 
to  separate  the  different  possible  effects,  subsidiary  experiments  must 
be  performed.  These  experiments  depend  on  the  forces  which  can  be 
applied  to  different  kinds  of  particles,, 

1.  Charged  Particles 

The  mass  spectrometer  should  remove  from  the  beam  any  charged  particle 
of  a  charge-to- mass  ratio  other  than  that  for  which  it  is  in  resonance,, 
However,  for  high  energy  ions  or  electrons,  the  transit  time  may  be  smaller 
than  that  required  for  tne  mass  filter  to  act.  For  ions,  this  energy  is 
greater  than  200  ev.  (Mass  spectra  have  been  taken  with  ion  energies  of 
up  to  400  ev  and  the  filter  is  effective  at  this  energy).  Thus,  for  ions  up  to 
200  ev  initial  energy,  the  mass  spectrometer  set  to  be  electronically  opaque 
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should  stop  all  ions.  Electrons  of  relatively  low  energy  will  have  sufficient 
velocity  to  traverse  the  spectrometer  without  being  filtered.  These  electrons 
would,  however,  not  be  able  to  strike  the  first  dynode  of  the  multiplier  unless 
they  had  at  least  3000  ev  of  energy.  Thus,  setting  the  spectrometer  to  be 
electronically  opaque  should  stop  all  charged  particles  of  less  than  200  ev. 

In  addition,  a  magnetic  field  of  approximately  500  gauss  can  be  applied 
externally  to  ths  spectrometer  housing  crosswise  to  the  path  of  the  particles. 
Such  a  magnetic  field  will  give  even  high  energy  ions  curvature  of  path 
following  the  equation: 


R  = 


144  mE 
eB  ^ 


1/2 


where  R  is  the  radius  of  curvature  of  the  trajectory  in  cm,  E  is  the  particle 
energy  in  ev,  B  is  the  field  in  gauss  and  e/m  is  the  charge  to  mass  ratio  in 
electronic  charge 'atomic  mass  units.  For  B  -  500  gauss,  E  ••  200  ev  and 
e/m  -  1/50  e/amu,  R  is  about  2  cm.  Since  this  500  gauss  field  could  be 
applied  over  an  area  encompassing  at  least  5  cm  of  particle  path  length,  it 
should  be  possible  to  stop  relatively  fast  ions  with  it  as  well  as  with  the 
mass  filter.  In  the  experiments  to  be  described  below  both  methods  are 
used. 


2.  Photons 


If  photons  are  emitted  from  the  point  of  impact  of  the  laser  light,  some  of 
them  will  strike  the  first  dynode  of  the  electron  multiplier  and  will  be  detected 
via  secondary  electron  effects.  The  photons  will  require  an  energy  greater 
than  the  work  function  of  the  beryllium -  copper  dynode.  This  energy  is  not 
well  known  but  should  be  equivalent  to  ultraviolet  radiation,  as  the  multiplier 
is  known  to  be  insensitive  to  the  white  hot  filament  in  the  ion  source  of  the 
mass  spectrometer.  Energy  analysis  of  the  photons  may  be  accomplished 
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by  passing  them  through  a  filter.  For  crude  analysis,  rather  simple  filters 

o 

may  be  used  such  as  quartz  (passes  wave-lengths  longer  than  2000A  ), lithium 

o  o 

fluoride  (cutoff  at  1100A  ),  barium  fluoride  (cutoff  at  1350A),  and  pyrex 

O 

(cutoff  at  2800A),  Experiments  are  discussed  below  using  these  filters  to 
give  the  approximate  spectrum  of  the  ultraviolet  radiation  emitted.  The 
spectrum  must  be  corrected  for  the  wavelength  dependence  of  the  photo¬ 
electric  yield  of  the  first  dynode  of  the  detector. 

3 .  Uncharged  Particles 

Uncharged  par.ides  of  atomic  dimensions  can  be  stopped  only  by  interposition 
of  a  physical  barrier.  Materials  such  as  quartz  may  be  used  which  will  pass 
ultraviolet  radiation  but  stop  uncharged  (as  well  as  charged)  species.  Since 
no  moderate  fields  will  deflect  uncharged  particles,  they  will  pass  through  the 
spectrometer  undeflected  and  will  be  detected  by  the  electron  multiplier  if 
the  kinetic  energy  is  sufficient  to  cause  secondary  emission  from  the  surface 
of  the  first  dynode  of  the  multiplier.  V/e  will  assume  that  this  will  be  true  for 
atoms  of  energy  greater  than  approximately  50  ev.  In  the  event  that  the  atoms 
are  excited,  however,  the  energy  of  excitation  will  be  available  for  secondary 
electron  ejection,  causing  detection  at  even  lower  kinetic  energies.  Since 
these  atoms  will  be  moving  at  considerably  slower  speed,  they  will  arrive 
quite  late  after  the  photon  pulse.  The  kinetic  energy  (not  including  energy  of 
excitation)  of  high-speed  atoms  can  be  inferred,  if  one  knows  the  mass  of 
the  atom,  from  the  time  of  flight  of  the  atom  from  its  source  to  the  detector. 
For  a  mass  of  m  amu  and  for  a  kinetic  energy  of  E  electron  volts,  the  time 
in  seconds  for  a  flight  of  d  cm  is: 

t  =  10'6  d  [  m/2E]  1/2 
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which,  for  m  “  1  amu,  d  -  60  cm  and  E  =  200  ev, results  in  t  -•  3  /usee . 

For  m  -  1  amu,  d  =  60  cm  and  E  =  50  ev,  t  =  6  /usee.  Thus,  the 
results  of  having  a  spread  of  ion  energies  from  0  to  200  ev  will  result  in 
a  peak  from  3  to  6  /u  seconds  after  the  initiation  of  the  pulse  of  high  speed 
atoms.  (No  consideration  of  atoms  slower  than  50  ev  is  taken  because 
the  secondary  electron  emission  is  very  small  for  energy  less  than  50  ev. ) 

It  is,  of  course,  impossible  to  accelerate  or  decelerate  the  atoms  so  they 
cannot  be  separated  by  electric  fields  as  are  the  ions  in  the  time -of- flight 
spectrometer. 

D.  STUDY  OF  THE  HTGH  SPEED  NEUTRAL  MOLECULES 

In  our  December  13th  progress  report  we  discussed  the  first  observation 
of  what  were  possibly  high-speed  neutral  molecules  emitted  from  the  laser 
target.  At  that  time  it  was  clear  that  there  was  some  sort  of  signal  which 
was  occurring  very  shortly  after  the  laser  burst, and  which  was  apparently 
unaffected  by  any  applied  fields. 

Since  that  time  considerable  work  was  done  to  determine  just  what  these  pulses 
were.  Our  conclusion  is  the  following:  When  a  Q- switched  laser  beam  strikes 
a  target  surface,  there  is  an  immediate  burst  of  ultraviolet  light  extending 

O 

down  to  wavelengths  of  the  order  of  11 00 A  accompanied  by  an  emission  of 
neutral  molecules  of  atomic  weight  of  1  or  2  amu  and  up.  These  neutral 
molecules  have  energies  of  the  order  of  100  ev. 

A  typical  trace  taken  from  the  data  is  shown  as  figure  1-1.  None  of  the  pulses 
shown  disappear  when  the  spectrometer  is  set  to  be  electronically  opaque  or 
when  a  crossed  magnetic  field  of  500  gauss  is  applied;  therefore,  one  may 
conclude  that  they  are  not  caused  by  charged  particles  of  either  sign.  These 
peaks  are,  therefore,  due  to  high  speed  atoms,  excited  atoms,  or  photons. 
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1, 


The  Photon  Pulse 


We  have  shown  that  the  first  peak  in  this  picture  is  a  photon  pulse.  The 
experiment  used  to  show  this  involves  the  use  of  filters  of  LiF,  BaF2» 
quartz,  and  pyrex.  Table  1-2  shows  the  height  of  this  peak  when  filters 
of  these  materials  are  interposed  between  the  target  and  the  electron 
multiplier  detector.  Note  that  this  is  not  meant  to  be  the  UV  emission  vs 
wavelength  but  the  product  of  the  UV  emission  times  the  quantum  efficiency 
of  the:  first  dynode  vs  wavelength.  It  should  be  noted  that  the  height  of  this 
peak  decreases  as  the  cutoff  wavelength  of  the  filter  is  made  longer.  The 
fact  that  there  is  some  amplitude  of  this  pulse  when  pyrex  glass  is  used  as 
a  filter  should  not  be  used  to  infer  that  there  is  radiation  of  wavelength 

O 

longer  than  2800A  .  The  filter  does  not  block  all  paths  from  the  target  to 
the  detector,  only  direct  ones;  so  there  may  we11  be  some  reflected  UV 
passing  around  the  filter.  The  data  show  that  the  height  of  the  first  pulse 
varies  monotonically  with  the  power  density  of  the  laser  illumination  of 
the  target.  The  fact  that  the  first  pulse  of  the  train  (1)  has  a  height  which 
varies  monotonically  with  the  illumination  density,  (2)  is  not  changed  by 
interposing  electric  or  magnetic  fields  between  the  target  and  the  detector, 

(3)  is  decreased  in  intensity  in  agreement  with  the  characterist.es  of  known 
UV  filters,  and  (4)  occurs  coincidentally  with  the  laser  burst,  establishes 
without  doubt  that  it  is  due  to  ultraviolet  light  emitted  by  the  target  upon 
irradiation  by  the  Q-switched  laser  beam. 

2.  The  High  Energy  Neutral  Molecule  Pulses 

Following  the  pulse  of  photons  in  figure  I  - 1 ,  are  three  other  pulses.  We 
have  shown  that  none  of  these  pulses  is  affected  by  (0  making  the  spectrometer 
electrically  opaque,  (2)  applying  a  crossed  magnetic  field,  or  (3)  applying 
positive  or  negative  potentials  to  the  target.  They  can  be  made  to  dissappeur, 
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TABLE  1-2 


TRANSMISSION  OF  THE  FIRST  PULSE  (PLOTON  PULSE) 

FOR  VARIOUS  FILTERS 


FILTER  MATERIAL 

APPROXIMATE 
UV  CUTOFF  (A) 

AMPLITUDE  OF  FIRST  PULSE 
AS  FRACTION  OF  AMPLITUDE 
WITH  NO  FILTER 

LiF 

1100 

0.7 

BaF  2 

1350 

0.5 

QUARTZ 

2000 

0.3 

PYREX 

2800 

.01-. 02 
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however,  by  interposing  any  of  the  UV  filters  described  above  between  the 
target  and  the  detector.  Furthermore,  they  are  delayed  in  time  from  the 
laser  burst  by  times  in  excess  of  2  Msec.  This  set  of  conditions  makes  it 
certain  that  these  pulses  are  due  to  high  energy  neutral  particles. 

The  first  of  these  pulses  can  be  shown  to  be  either  due  to  neutral  hydrogen 
atoms  with  kinetic  energies  of  the  order  of  100  ev,  or  some  other  neutral 
particle  with  even  higher  energy.  As  we  have  shown  in  Section  I-C-3  the 
time  of  flight  of  a  neutral  particle  of  mass  m  amu  and  with  energy  E  ev  over 
a  distance  d  cm  is,  in  seconds: 

t  =  10"6  d  Em/2E] 

Thus,  for  a  given  time  of  flight,  an  increase  in  m  must  be  accompanied  by  a 
proportional  increase  in  E.  Since  the  first  of  the  pulses  we  have  called  the 
neutral  particle  pulses  negins  to  arrive  at  the  detector  only  5  Msec  after  the 
laser  burst,  and  since  the  detector  is  60  cm  from  the  target,  the  particle 
must  have  either  a  mass  of  1  amu  (the  smallest  mass  possible  for  a  neutral 
particle)  and  a  kinetic  energy  of  70  ev,  or  a  mass  greater  than  1  and  even 
larger  energy.  The  data  reported  does  not  permit  any  other  possibilities. 

3.  M  ass  Discrimination  Efforts  in  the  High  Seeed  Neutral  Molecule  Pulses 


it  is  not  possible  to  accelerate  atoms  by  fields.  As  a  result,  it  is  net  possible 
to  determine  separately  the  energy  and  mass  of  the  high  speed  atoms  as  can 
no  done  with  the  high  speed  ions.  It  is  possible,  however,  to  determine  a 
minimum  mass -maximum  energy  crossover  for  a  given  peak  and  to  match 
peaks  from  a  given  laser  burst  for  energy  and  mass  possibilities. 

This  sor  t  of  analysts  is  carried  out  in  Sec  tion  III  using  typical  da* a  taken 
with  this  apparatus. 
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STUDIES  ON  HIGH  ENERGY  MOLECULES 


IT 


High  energy  molecules  have  been  observed,  thus  far,  emanating  from  tungsten, 
single  crystal  nickel,  and  titanium.  Tungsten  was  studied  because  it  was  the 
metal  on  which  the  bulk  of  our  work  was  performed.  When  the  high  energy 
atoms  were  detected  emanating  from  tungsten,  it  became  clear  that  some 
other  metal  would  have  to  be  used  to  determine  the  species  since  it  is  well 
known  that  all  kinds  of  gas  can  reside  on  a  surface  of  polycrystalline  tungsten. 
It  was  decided  to  use  as  our  target  two  kinds  of  surfaces  -  single  crystal 
(100)  nickel  and  polycrystalline  titanium. 

Work  by  D.  Lichtman  in  our  laboratory  has  shown  Uiat  single  crystal  (100) 
nickel  seems  to  have  hydrogen  as  its  only  gas  coverage.  When  different  gases 
are  admitted  into  his  vacuum  system  and  electron  induced  ion  desorption  is 
attempted,  only  ions  resulting  from  hydrogen  desorption  have  ever  been 
observed.  Thus,  it  was  hoped,  only  one  species  of  high  speed  atom  or 
molecule  would  be  detected  in  the  laser  induced  desorption  experiment  with 
(100)  nickel,  and  that  species  would  be  of  mass  1  or  2  amu.  With  only  one 
species  of  desorption,  the  analysis  of  the  energy  of  the  atom  or  molecule 
would  be  quite  simple.  Titanium,  on  the  other  hand,  can  be  covered  by  a 
wide  range  of  gases.  In  fact,  the  surface  is  so  active  chemically  that  all 
but  the  rare  gases  are  chemically  bonded  to  the  surface.  Of  these  gases, 
only  hydrogen  is  bonded  reversibly.  It  was  felt  that  with  titanium  as  a 
target,  the  usual  array  of  high  speed  atom  pulses  would  be  observed,  but 
that  only  one  of  them  -  the  first  in  time  -  could  be  made  to  vanish  as  the 
target  was  heated  to  a  temperature  in  excess  of  30CF  C. 
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The  studies  on  the  nickel  target  were  inconclusive.  When  the  laser  beam 
was  allowed  to  strike  the  target,  peaks  were  observed  which  could  not  be 
interpreted  in  terms  of  only  one  species  of  gas  on  the  target.  When  gas 
desorption  data  were  taken,  it  was  found  that  both  hydrogen  gas  and  carbon 
monoxide  gas  were  desorbed.  When  the  surface  was  heated  to  as  high  as 
70(F  C  (a  temperature  above  that  at  which  hydrogen  was  found  by  Lichtman 
to  desorb  from  the  surface),  the  first  peak  of  the  high-speed  atom  gi  onp 
vanished  at  times  and  remained  at  others.  The  crystal  was  removed  from 
the  vacuum  system  and  subjected  to  X-ray  diffraction  analysis.  It  was 
found  that  the  region  struck  by  the  laser  had  been  melted  and  refrozen  in  a 
polycrystalline  form.  The  data  were  therefore  totally  in  agreement  with 
what  could  have  been  the  case  T  r  polycrystalline  material. 

When  titanium  was  used  as  a  target,  the  results  were  somewhat  more 

conclusive.  The  data  had  generally  the  appearance  of  3  or  4  peaks.  The 

first  peak  did  in  fact  come  and  go  as  a  function  of  the  temperature  of  the 

target,  and  the  other  peaks  did  not.  Figure  1-2  shows  typical  data  for  a 

titanium  target.  In  both  I-2a  and  I-2b,  time  goes  to  the  right.  In  each 

case  the  top  trace  (signal  extending  downward)  has  a  speed  of  2  ti sec/cm 

_  3 

and  a  sensitivity  of  10  amps/cm.  The  lower  trace  (signal  extending  up- 

_  3 

ward)  has  a  speed  of  5  /jsec/cm  and  a  sensitivity  of  4  x  10  amps/cm. 
Figure  I-2a  is  for  a  cold  titanium  target.  Figure  I-2b  is  for  a  hot  titanium 
target.  The  only  important  difference  between  them  is  a  disappearance  of 
the  second  peak  from  the  left  in  figure  I-2b. 

It  was  not  possible,  however,  to  establish  the  presence  or  absence  of 
gaseous  hydrogen  desorption  as  a  function  of  target  temperature  wiiich 
correiated  with  the  high  speed  atom  data.  (The  reflected  laser  energy 
for  these  runs  was  inadvertently  allowed  to  strike  an  interior  metal 
surface  instead  of  the  exit  window  as  had  been  designed.  Gas  from  the 
surface  struck  by  this  unfocused  but  still  intense  beam  could  easily  have 
confused  the  data.  Work  is  now  underway  to  correct  this.  )  With  the 
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(a)  Cold  Target 


(b)  Hot  Target 


Figure  1-2  -  Photon  Pulse  and  High  Speed  Molecule  Pulse  from  a  Titanium 
Target.  In  each  photograph  top  trace,  negative  signal  downward;  bottom 
trace,  negative  signal  upward.  Time  progresses  to  the  right.  Upper  trace: 

2  microseconds /cm.  Lower  trace:  5  microseconds /cm.  (a)  Cold  Target 
(b)  Hot  Target. 
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titanium  data  it  is  now  clear  that  the  presence  and  absence  of  the  first 
peak  of  the  high-speed  atom  group  can  be  correlated  with  the  temperature 
of  the  titanium  target.  Although  the  actual  temperature  at  which  this 
peak  disappears  is  not  yet  known,  it  is  in  a  range  which  correlates  with  the 
degassing  of  the  titanium  hydride.  The  other  peaks  remain  unaffected  by  the 
thermal  treatment  of  the  target.  This  is  in  agreement  with  the  interpretation 
of  these  peaks  as  due  to  other  nondesorbing  gases,  and  even  to  the  base 
material  itself. 

F.  THE  TOTAL  NUMBER  OF  MOLECULES  EMITTED  FROM  THE  TARGET 

Ar  estimate  may  be  made  of  the  total  number  of  atoms  of  a  particular  species 

emitted  from  the  target  during  a  laser  burst  if  one  knows  the  value  of  Y  for 

the  species  on  beryllium-copper  and  the  spatial  distribution  of  the  high  speed 

atoms,  in  our  case  neither  of  these  is  known  with  even  crude  precision. 

However,  we  ccn  say  that  for  the  energy  range  of  interest,  Y  is  in  excess  of 

0.  01  and  that  the  spatial  distribution  is  approximately  isotropic.  Using 

these  two  factors,  and  the  height  and  time  duration  of  the  pulse  coupled  with 

4 

the  gain  of  the  multiplier  (of  the  order  of  10  ),  we  can  arrive  at  the  total 

number  of  atoms  emitted.  A  typical  pulse, pulse  1  of  figure  1-1,  is  of  the 

order  of  0.  01  volt  high  with  an  input  resistor  of  500  ohms.  Simple  calculation 

6 

shows  that  this  signal  is  equivalent  to  10  electrons  emitted  from  the  first 

dynode  of  the  electron  multiplier.  With  a  value  of  Y  equal  to  0.  01,  this 
8 

implies  that  10  high  speed  atoms  arrived  at  the  detector.  The  detector  has 

2 

an  active  area  of  less  than  1  cm  and  is  60  cm  from  the  target,  If  one 
assumes  that  the  emission  of  high  speed  neutrals  is  isotropic  over  a  hemi¬ 
sphere,  the  total  number  of  nigh  speed  atoms  emitted  is  3.  6  x  10**.  This 
number  is  quite  high  when  compared  with  the  number  of  high  speed  ions 
observed.  It  can  be  reduced  by  two  orders  of  magnitude  by  changing  the 
estimate  of  gamma  from  0.  01  to  l,  a  not  unreasonable  difference.  It  can 
further  be  reduced  by  assuming  a  non  -  isotropic  distribution  of  atoms. 

(What  shape  distribution  should  be  used,  however,  is  dependent  on  the 
formation  mechanism  of  the  high  speed  molecules.  ) 
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II.  tom  emission 


The  work  on  ion  emission  for  this  period  has  been  concentrated  primarily  on 
the  interpretation  of  the  pulse  shapes.  On  the  experimental  side  we  have 
obtained  photographs  of  pulses  with  enough  time  resolution  to  be  useful  in  the 
analysis.  We  have  also  measured  the  dependence  of  the  time  of  flight  of 
several  mass  peaks  on  the  accelerating  potential  used  in  the  TOF  spectrometer*. 
These  results  and  their  use  in  finding  the  velocity  distribution  of  the  ions  are 
presented  in  this  section.  The  interpretation  of  data  is  presented  in  Section  III. 

The  ion  emission  from  a  thin  platinum  target  was  also  measured.  The  results 
seem  to  indicate  that  the  mechanism  of  ion  acceleration  involves  a  direct 
interaction  of  the  emitted  gas  with  the  laser  beam.  The  emission  itself  is 
qualitatively  similar  to  that  of  tungsten. 

A.  ION  EMISSION  FROM  PLATINUM 

The  purpose  of  measuring  ion  emission  from  platinum  is  twofold:  First,  it 
provides  us  with  a  different  surface  to  compare  with  tungsten;  second,  it 
allows  us  to  check  on  whether  or  not  the  laser  beam  is  directly  involved  in 
the  foi'mation  of  the  ions  observed. 

Figure  II- 1  shows  a  typical  spectrum  obtained  from  platinum  when  irradiated 
on  the  surface  facing  the  spectrometer.  The  spectrum  is  similar  to  that 
obtained  from  tungsten  surfaces  in  that  the  major  peaks  are  still  Na  and  K+. 
Other  peaks  are  present  but  they  are  much  smaller.  On  the  other  hand,  while 
K+  is  the  largest  peak  in  the  emission  from  tungsten,  the  Na+  peak  is  largest 
for  platinum.  The  ion  energies  from  both  surfaces  are  similar  for  equivalent 
laser  outputs.  In  general  the  emission  from  the  two  surfaces  (platinum  and 
tungsten)  is  quite  similar  except  for  the  relative  abundances  of  different 
s  pec  ies . 
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Figure  II- 1  -  Typical  Mass  Spectrum  from  Platinum  Target  in  TOF 
Spectrometer.  Time  increases  from  right  to  left.  Lower  trace  starts  with 
laser  pulse  (2  us/cm,  .  1  v/cm).  Upper  trace  is  an  expanded  portion  cf  the 
lower  trace  triggered  by  the  mass  23  pulse  (.  5 u  sec/cm,  .  05  v/cm. )  The 
spectrometer  potentials  are:  V  =  +1000v.  ;  v'  =  +1022v. 
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The  platinum  target  used  in  these  experiments  was  5  mit  rons  thick. 

Calculation  of  the  bulk  temperature  upon  irradiation  sh^ws  that  the  temper¬ 
ature  gradient  is  negligible  across  the  thickness  of  the  target.  Hence,  if 
the  mechanism  for  ion  production  and  acceleration  did  not  involve  a  direct 
interaction  with  the  laser  beam,  one  should  obtain  essentially  the  same  ion 
emission  regardless  of  which  side  of  the  target  was  irradiated. 

The  results  obtained  from  irradiation  of  the  back  side  of  the  target  were  all 
negative.  The  absence  of  ions  was  observed  in  experiments  that  consisted  of 
alternating  the  laser  beam  between  the  front  and  back  surfaces,  with  emission 
observed  always  from  the  front  surface.  These  results  are  consistent  with  the 
hypothesis  that  the  creation  and  acceleration  of  ions  is  associated  with  a 
direct  interaction  with  the  laser  beam  and  is  inconsistent  with  the  hypothesis 
that  ions  are  produced  by  a  shock  wave  traversing  the  target. 

B.  ELECTRIC  FIELD  DEPENDENCE  OF  TIME  OF  FLIGHT 

The  time  of  flight  was  measured  as  a  function  of  field  for  the  major  ion 
peak  of  tungsten  and  platinum  surfaces.  This  information  is  important  in 
determining  the  validity  of  mass  assignments  which  are  based  on  calculations 
that  assume  independent  particle  behavior  of  the  emission  from  the  metal.  In 
the  analysis  of  pulse  shapes  to  arrive  at  a  velocity  distribution  (which  is 
discussed  later  in  this  report),  it  is  assumed  that  the  time  of  flight  can  be 
written  as 


t 


vd(w) 


(II- 1) 


where  i 


effective  length  of  drift  region,  and 

vdM  -V5  k»-v,o)]  +w2 


( II- 2 ) 
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is  the  velocity  in  the  drift  tube  of  singly  charged  ions  of  mass  m  with  initial 
velocity  w  and  accelerated  by  a  potential  difference  V(xj)-V(o).  The  above 
express!'  rnvolve  two  assumptions:  first,  that  the  time  spent  outside  of  the 
drift  tube  is  small  compared  to  the  drift  time;  second,  that  the  particles 
emitted  by  the  target  behave  as  independent  particles  in  their  reaction  to 
externally  applied  fields. 

The  last  assumption  is  particularly  crucial  not  only  for  the  analysis  of 
distributions,  but  also  to  test  the  validity  of  our  conclusions  regarding  ion 
masses.  We  have  therefore  measured  the  time  of  flight  as  a  function  of 
applied  field.  If  the  above  assumptions  are  justified,  we  can  combine 
Equations  (II- 1)  and  (II- 2 )  to  obtain 


t-2  ,  2£_ 

o 

V(xd)  -  V(o) 

+  It) 

mi. 

_ 

Figures  II- 2  and  II- 3  chow  the  measured  arrival  times  as  a  function  of 
applied  voltage  for  the  largest  mass  peaks  obtained  from  tungsten  and 
platinum;  they  correspond  to  masses  39  and  23  respectively.  The  points 
represent  the  measured  values  and  the  curves  represent  the  best  fit  to 
the  data  in  the  least  squares  sense. 

It  can  be  seen  that  the  data  fits  the  straight  line  predicted  by  Equation  II- 3 
quite  well.  We  take  this  to  be  evidence  of  the  fact  that  the  plasma  consisting 
of  the  particles  emitted  bv  the  target  is  very  short  lived;  i.e.,  macroscopic 
charge  neutrality  is  not  maintained  over  distances  comparable  to  the  dimen¬ 
sions  of  the  accelerating  gap.  This  measurement,  therefore,  is  another 
piece  of  evidence  which  justifies  the  assignment  of  ion  masses  on  the  basis 
of  the  times  of  flight  measured  in  the  spectrometer. 
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Figure  II-  2  -  Time  of  Flight  vs.  Aeeelerating  Potential  for  the  Na  (22) 
Peak  from  pl;*Mnum  Target 
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One  further  question  that  must  be  answered  before  proceeding  with  the 
analysis  of  the  measured  pulse  shapes  is:  What  is  the  effective  length  l 
that  should  be  used  in  calculating  t  in  Equation  (II- 1)?  Neglecting  initial 
ion  energies  as  compared  to  those  given  by  the  fields,  it  is  easy  to  show 
that  the  effective  "drift  length"  of  the  spectrometer  is  ~  105  cm.  The 
argument  is  as  follows:  Let  x  =  a  in  the  accelerating  gap  and  v^w)  be  the 
velocity  at  the  end  of  the  gap  (entrance  to  drift  tube);  w  is  the  initial  particle 
velocity.  Then, 

v(w)  =  at  +  w . 


If  we  let  t^  be  the  time  required  to  cross  the  gap. 


v^(w)  -  at^  +  w 


a 

v  (w)  =  — —  I  v(w)dt  =  j  a  + 
ld  '  o 


w 


and 


xd  '  2-atd2+wtd  =  0  <w,td 


(11-4) 


We  assume  that  at  ^  >>  w;  i.e.,  the  energies  supplied  by  external  fields  are 
much  larger  than  the  initial  energies.  Then 


v 


d 


(II"  5) 


Using  (II- 5)  in  (II-  4),  we  obtain 


=  Vd 


(Il-fi) 


Hence,  since  we  wish  to  treat  the  spectrometer  as  if  the  particles  had  a 
constant  velocity  v^  while  traversing  it,  (II- 6)  gives  the  prescription  for 
handling  the  regions  where  fields  are  applied.  For  the  particular  field 
configurations  used  in  the  spectrometer^,  equation  (II- 6)  says  that  we  must 
double  the  actual  length  of  the  region.  Using  this  prescription  we  obtain  an 
effective  drift  length  X — -105  cm. 

With  the  drift  length  we  can  now  check  the  mass  assignments  of  the  peaks  of 

Figures  II- 2  and  11-3.  The  slopes  corresponding  to  the  least  squares  fit  are 

+6  +5 

4.  7  x  10  0  for  the  peak  from  tungsten  and  7.  8  x  10  '  for  the  peak  from 

platinum.  The  slopes  calculated  for  masses  39  and  23  usii.y  Equation  (II - 3 ) 
4*6  +g 

are  4.  65  x  10  and  7.  89  x  10  respectively.  The  agreement  between  the 
experimental  and  calculated  vaiues  is  another  confirmation  of  the  correctness 
of  mass  assignments. 


C.  EXTENSION  OF  PULSE  SHAPE  ANALYSIS 

In  the  last  semi-annual  report,  *  it  was  shown  that  for  a  monoenergetic 
pulse  of  ions  of  a  given  shape,  the  detector  signal  in  the  time-of-flight 
spectrometer  would  reproduce  the  shape  of  the  pulse.  In  this  report  we 
extend  the  proof  to  any  parallel-plate  structure  and  to  an  arbitrary  pulse 
shape.  The  one- dimensional  equation  of  continuity 


+  1  ~J 

31  x  v  d  t 


0 


governs  the  flow  of  charge  in  Mu*  spectrometer. 


Its  La  place  transform 


.  a 
is 


lLl 

(|v 


+  _ 


V 


0 


(II  7) 


Equation  (II- 7)  has  a  solution 


j  (x,  s)  =  j(o,  s)  exp  I  -s  I  --  =  L  [  J  (x,  t)  ] 


J  v 


where  j(o,  s)  L  [  J(o,  t)  1, 


(11-0) 
(II- 9) 


i.e.,  j(o,  s)  is  the  Laplace  transform  of  the  boundary  conditions.  In  general, 
v  may  or  may  not  be  constant,  depending  on  whether  or  not  external  fields  are 
present.  However,  it  is  a  well  known  property  of  Laplace  transforms  that  if 
(II- 9)  is  true,  then"* 


j  (o,  s )  exp 


1 

72 

lx 

_ J 

r 

J(o,t  -  f*L)' 

J  V  j 

J  V 

(II- 10) 


Comparing  (II - 8 )  and  * II  ■  10),  we  find  that  for  an  arbitrary  parallel-electrode 

V 

geometry,  or  for  that  matter  any  structure  for  which  —  -  0,  the  convection 

at 

current  at  any  point  is  simply  the  current  at  the  target  delayed  by  the  time  of 
flight  from  target  to  detector.  Notice  that  the  proof  is  independent  of  any 
assumptions  about  pulse  shapes. 


The  above  proof  suggests  a  possible  new  structure  that  can  be  used  to 
obtain  information  on  both  the  density  o(o,  w,  t)  at  the  target  and  the  initial 
ion  velocity  w.  Consider  a  triode  consisting  of  a  target,  an  accelerating 
grid  followed  by  a  drift  region,  and  a  collector.  Using  the  results  of  the  last 
semi-annual  report*,  one  finds  that  if  the  detector  is  a  plain  electrode,  then 
the  detector  current  I  satisfies 


—  s  —  J(o,i)  (n-  ii) 

dt  d 

where  w  :  initial  ion  velocity 
A  r  detector  area 

d  -  separation  between  accelerating  grid  and  detector. 
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On  the  other  hand,  if  one  uses  an  electron  multiplier  detector,  the  current  is 
given  by 


I(t)  =  gA  J(o,  t);  g  =  multiplier  gain.  (11-12) 

The  above  arrangement  can  be  used  as  follows:  If  a  target  emitting  particles 
of  only  one  mass  is  used,  the  ratio  of  the  two  signals  given  by  (II- 11)  and  (11-12) 
gives  the  initial  ion  velocity  w  as  a  function  of  time.  The  two  signals  can  be 
derived  from  different  electrodes  in  a  two  stage  multiplier  detector.  Knowledge 
of  w  allows  one  to  compute  the  ion  density.  In  practice  the  emission  is  not 
moncenergetic,  and  it  would  be  necessary  to  bias  the  detector  to  admit  only  a 
small  range  of  velocities.  The  single  mass  requirement  may  possibly  be 
satisfied  either  by  using  laser  powers  that  are  high  enough  to  evaporate  a 
considerable  amount  of  the  target  material  per  shot,  or  by  using  a  (100)  face 
of  single  crystal  nickel  that  can  be  covered  with  hydrogen  exclusively. 


III.  INTERPRETATION  OF  DATA 


This  portion  of  the  report  covers  the  interpretati  3n  of  the  experimental  results 
described  in  the  previous  two  sections*  It  includes  the  work  that  is  basically 
theoretical  or  computational.  There  are  three  main  topics  considered  in  this 
section:  The  first  involves  processing  the  measured  pulse  shapes  of  emitted 
particles  to  infer  their  original  velocity  spectrum.  The  second  involves 
extensions  of  work  described  in  the  semi-annual  report  ;  namely,  the  cal¬ 
culation  of  the  heating  effects  that  arise  from  absorption  of  the  laser  light 
by  free-free  transitions  in  the  blow  off  material.  The  third  is  a  discussion  of 
the  identification  of  the  masses  and  energies  of  the  neutral  particles  described 
in  Section  I.  While  the  identification  is  not  definite,  some  tentative  hypotheses 
are  made. 

A.  DETERMINATION  OF  VELOCITY  DISTRIBUTION 

The  calculation  of  the  ion  velocity  distribution  from  the  measured  pulse 
shapes  is  based  on  the  analysis  given  in  Appendix  I.  We  will  review  briefly 
here  those  parts  that  are  rel°vant  to  this  discussion  and  indicate  how  the 
data  is  processed  Basically  we  assume  that  the  ion  pulse  from  the  target 
is  a  6-function  in  time,  i.  e. ,  we  assume  that  the  finite  width  of  the  pulse  at 
the  detector  is  due  to  a  distribution  in  initial  ion  velocities.  This 
assumption  is  justified  on  the  basis  that  the  width  of  the  ion  pulse  at  the 
detector  is  much  longer  than  that  of  the  laser  pulse  during  which  one  would 
expect  the  ion  emission  to  occur.  The  above  assumption  establishes  a 
one-to  one  correspondence  between  time  of  arrival  and. initialion  velocity  with  the 
time  of  arrival  of  zero  velocity  ions  corresponding  to  the  end  of  the  detector 
pulse. 

Now,  the  velocity  v^(w)  of  a  singly  cha  ged  ion  of  mass  m  with  initial 
velocity  w,  accelerated  through  a  potential  difference  V(x^)-V(o),  is 

vrf(w)  =  "W— ^Tv(xd)  ~  V(o)j+ w4"  (III- 1 ) 
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1 


Furthermore,  since  most  of  the  time  of  flight  t  is  spent  in  the  drift  tube, 
we  approximate  it  by 


t 


l 


vd(w) 


(III- 2 ) 


where  i  is  an  effective  drift  length  that  corrects  for  time  spent  outside  the 
drift  tube.  The  justification  for  this  approximation  was  given  in  Section  II- B. 
Notice  that  (III- 2 )  implies  the  6-function  assumption  for  the  emission  pulse. 
From  (III- 1 )  and  (III- 2 )  we  get 


dt  _  Xw 


dw 


vd3(w) 


(III  -  3 ) 


Now,  if  we  let  N(w)  represent  the  distribution  of  velocities  of  the  ions  at  the 
target,  conservation  of  charge  requires  that 


N(w)  dw  = 


J(t)  dt 


(III- 4) 


where  J(t)  is  the  signal  from  the  detector;  k  is  the  constant  transmission 
coefficient  of  the  spectrometer.  We  can  rewrite  (III- 4 )  as 


00 


O 


kN(w)  -  J(t)  — 
dw 


dw  =  0 


from  which  we  find  the  prescription  for  calculating  N(w), 

N(w)  =  k~  1  J(t)  —  . 

dw 


(III- 5 ) 
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The  application  of  the  above  to  the  data  is  as  follows:  Values  of  J(t) 

measured  at  equal  intervals  are  stored  in  a  computer.  The  end  of  the  pulse 

is  taken  to  be  the  time  of  arrival  of  the  ions  with  zero  initial  velocity  (w  =  0). 

A  constant  initial  velocity  increment  dw  is  chosen  and  added  to  w.  With  the 

new  w  used  in  (III-  1),  (III- 2 )  and  (III- 3),  new  value's  of  t’aridf  tft/dVr  afie'caletilslted. 

These  values  are  then  used  to  calculate  N(w)  using  (III- 5 )  and  values  of  J(t) 

interpolated  from  the  data  stored  in  the  computer.  The  process  continues 

2 

with  a  new  increment  dw.  At  each  step  the  quantities  wN(w)  and  w  N(w)  are 
calculated  and  stored  cumulatively  to  calculate  the  average  velocity  w  of  the 
center  of  mass,  as  well  as  the  random  energy  about  the  center  of  mass.  The 
average  velocity  is  calculated  from  Equation  (19a)  of  Appendix  I: 

w  =  MMj!  (Ill- 6) 

N 

where  N  =  ?  N(w)  dw.  The  random  energy  of  Equation  (19b)  of  the  Appendix  is 
calculated  using  the  following  simplification: 

—  T  (w-w)2  N(w)  dw  =  —  T  (w2-2w  w  +  (w  )2)  N(w)  dw  =  —  Ew2N(w)  dw  -  (w  )2 
N  N  N 

Using  the  above  method  one  eliminates  the  need  to  calculate  w  first,  and  then 

-  2 

repeat  the  calculation  to  form  (w  -  w)  at  each  stage. 

A  sample  of  the  computer  output  of  a  typical  calculation,  corresponding  to  the 
results  shown  in  Figure  III- 3  is  given  in  Table  III—  1. 

A  word  of  caution  seems  in  order  at  this  point.  While  the  calculations 
described  above  are  straightforward,  it  is  important  to  choose  a  small 
value  of  the  velocity  increment  dw.  If  dw  is  too  large,  the  measured  pulse 
shape  J(t)  changes  drastically  during  the  corresponding  time  increment  dt, 
and  the  resulting  velocity  distribution  appears  deceivingly  symmetrical  In 
fact,  in  some  of  our  earlier  calculations  of  velocity  distributions,  the 
results  looked  like  a  good  fit  to  a  Maxwellian  because  of  this  effect. 
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TABLE  III-  1 

TYPICAL  CALCULATION  OF  N(w)  FROM  MEASURED  PULSE  SHAPE  J(t) 


w 

T 

J  (T) 

N  (W) 

*  0.000L-99 

4  1.00E+05 

*  8.2081E-C6 

•S.2078E-06 

41.6155E-03 

48.14901-18 

+2.00E+05 

+8.2071E-06 

46 .461  IE-03 

4  6.5165E-17 

+3.00E+05 

+8.2058E-06 

41.4534E-02 

+2.1978E-16 

+4.00E+05 

+8.2041E-06 

42.5830E-0 2 

+5.2046E-16 

+5.00E+05 

+8.2018E-06 

44.0343E-02 

+1.0152E-15 

+6.00E+05 

*  8. 1990E-06 

*  5.6064E-02 

+1.7517E-15 

+  7.00E+05 

+8. 1958E-06 

47.8985F-02 

+2. 7766E-15 

+8.00E+05 

+8. 1920E-06 

4 1 .0309E-01 

44 . 1 362E-15 

J  9.00E+05 

+8.1877E-06 

41.3037E-01 

+5. 8755E-15 

+1.00E+06 

+8. 1830E-06 

-41.6081E-01 

+8.0386E-15 

+1.10E+06 

+8. 1777E-06 

+1.9440E-01 

+1.0668E-14 

+1.20E+06 

+8. 1720E-06 

42.3111E-01 

+  1.3807E-14 

+1.30E+06 

+8. 1658E-06 

+2.7092E-01 

+1.7494E-14 

+1.40E+06 

+8. 1591E-06 

+3. 1382E-01 

+2.1769E-14 

+1.50E+06 

+8. 1519E-06 

+3.5838E-01 

+2.6566E-14 

+1.60E+06 

+8. 1443E-06 

+4.0037E-01 

+  3.1568E-14 

+1.70E+06 

+8.1361E-06 

44.4495E-01 

+3.7164E-14 

+1.80E+06 

+8.1276E-06 

+4 .9207E-01 

+4.3381E-14 

+1.90E+06 

+8.1185E-06 

+5.4 1 74E-01 

+5.0244E-14 

+2.00E+06 

+8. 1090E-06 

+5.9390E-01 

+5.7778E-14 

+2.10E+06 

+8.0991E-06 

+6.4855E-01 

+6.6005E-14 

+2.20E+06 

+8.0887E-06 

+6.8709E-01 

47.2976E-14 

+2.30E+06 

+8.0778E-06 

+7.2677E-01 

+8.0375E-14 

+2.40E+06 

+8.0665E-06 

+7.6804E-01 

+8.8262E-14 

+2.50E+06 

+8.0548E-06 

+8  1088E-01 

+9.6646E-14 

+2.60E+06 

+8.0427E-06 

+8.6188E-01 

+1.0635E-13 

+2.70E+06 

+8.0302E-06 

49.6520E-01 

+1.2310E-13 

+2.80E+06 

+8.0172E-06 

+1.0719E+00 

41.4109E-13 

+2.90E+06 

-48.00  38E-06 

+1.1819E+00 

+1.6033E-13 

+3.00E+06 

+7.9901E-06 

+1.2953E+00 

+1.3082E-13 

+3.10E+06 

+7.9759E-06 

+1.7226E+00 

42.4718E-13 

+3.20E+06 

+7.96I4E-06 

+  2. 1750E+00 

+3.2040E-13 

+3.30E+06 

+7.9465E-06 

+2.6391E+00 

+  3.9867F.-13 

+3.40E+06 

47.9312E-06 

+  3. 1553E+00 

44 . 8826E-13 

+3.50E+06 

+7.9155E-06 

+3.8143E+00 

+6.0400E-13 

+3.60E+06 

+7.8995E-06 

44.4882E+00 

+7.2660E-13 

43.70E+06 

+7.8831E-06 

+5.1769E+00 

+8.5603E-13 

+3.80E+06 

+7.8664E-06 

+6.6116E+00 

+1.1156E-12 

+  3.90E+06 

+7.8494E-06 

4  8.2332E+00 

4l  .4166E-12 

44.00E+06 

+  7.8320F.-06 

+9 .8860E+00 

+1.7330E-12 

+4.10E+06 

47.8143E-06 

4 1 .0335E+01 

+1.8446E-12 

+4.20E+06 

+  7. 7963E-06 

41.0071E401 

41.8287E-12 

44.30E406 

t-7. 7780E-06 

49.8038E+00 

+1.8096E-12 

+4.40E+06 

+7 . 7594E-06 

48.8153E+00 

+1.6530E-12 

+4.50E+06 

4  7. 7405E-06 

47.3629E400 

41.4018E-12 

44.60E406 

4  7. 7213E-06 

45.8888E+00 

41 . 1375E-12 

44 . 70E+06 

47.7019E-0'' 

44.0512E+CO 

47.9358E-13 

44.80E406 

+7.6822E-06 

42.0663E400 

44 . 1020E-1 3 

44.90E406 

47.6622E-06 

'  5 . 5237E  -02 

4 1 . 1 107E-14 

45.00E4Q6 

4  7.6420E-06 

t0.0000r-99 

t  0  .0000E  -99 

N 

W 

(W-W)? 


t  2.024  7  7652E  -06 
t  3.9  731595  31  *06 
*  3. 543040601. *-11 
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Ion  Velocity  Distribution 


Several  of  the  expanded  photographs  of  mass  peaks  obtained  in  the  TOF 

spectrometer  have  been  subjected  to  the  analysis  described  above  in  an 

attempt  to  find  the  velocity  distribution  of  the  ions.  Before  analyzing  the 

data,  we  determined  experimentally  that  the  time  of  flight  can  be 

reasonably  approximated  by  an  expression  such  as  Equation  (III-2),  and 

determined  an  effective  length  to  be  used  in  that  equation.  Those  results  are 

reported  in  Section  II.  Figure  III-  1  shows  mass  pulses  obtained  from  a 

tungsten  target.  From  left  to  right  they  correspond  to  masses  23  and  39 

respectively.  The  accelerating  potential  applied  to  the  target  was  V  = 

+2000  volts.  The  velocity  distribution  obtained  from  the  analysis  of  the  mass 

pulses  are  given  in  Figures  III- 2  and  III  3  for  the  mass  39  and  mass  23 

peaks  respectively.  The  mass  39  pulse  has  a  center  of  mass  velocity,  w 

3.2  x  10  cm /see,  corresponding  to  a  directed  energy  of  the  center  of 

mass  of  207  ev  .  The  random  energy  assuming  equipartition,  is  16  ev. 

_  £» 

The  corresponding  parameters  for  the  mass  23  pulse  are  w  r  4.  0  x  10 
cm/sec,directed  center  of  mass  energy  of  190  ev,  and  random  energy  of  13 
ev .  In  both  Figures  III-- 2  and  III- 3  we  have  plotted  the  Maxwellian  distribution 

-  -  rn  (w  -  w)2 
L  2kT 

with  N,  w  and  T  taken  from  the  calculated  distributions. 

Qualitatively  it  appears  that  the  calculated  velocity  distributions  of  Figures 
III  -  2  and  III- 3  do  net  fit  a  Maxwellian  well,  for  they  are  asymmetric  and 
obviously  peaked  at  high  velocities.  Also,  the  maximum  energies  of  the  ions 
in  the  calculated  distributions  are  312  ev  for  the  mass  39  pulse  and  288  ev 
for  the  mass  23  pulse.  These  energies  are  about  100  ev  higher  than  the 
measured  ion  energies. 


N  (w)  =;  N 
m 


m 


2  nk  T 


1/  z 


exp 
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Figure  III- 1  -  Na+  and  k+  Pulses  from  Tungsten  Target  with  an  Accelerating 
Potential  of  +2000v.  Time  increases  to  left.  Lower  trace  triggers  with  laser 
pulse  and  sweeps  at  2  usec/cm.  Upper  trace  triggers  on  Na+  peak  and  sweeps 
at  0.  5  us /cm. 
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0 0  *  io‘6(cm/sec) 

Future  III-2  -  K+3g  Velocity  Distribution  N3g(w)  Calculated  from  Data  of 
Figure  III-  1 .  Maxwellian  Fit  Included  for  Comparison. 
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ARBITRARY  UNITS 


20 


Figure  III- 3  -  Na^.j  Velocity  Distribution  N^.^Cw)  Calculated  from  Data  of 
Figure  III- 1 .  Maxwellian  Fi*  Included  for  Comparison. 
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A  more  careful  examination  of  the  facts*  however,  suggests  that  the 
agreement  may  be  better  than  it  appears:  First,  the  distributions 
correctly  predict  that  the  ions  in  the  mass  39  pulse  are  more  energetic  than 
those  of  the  mass  23  pulse,  and  predicts  a  difference  of  —  30  ev  for  the 
maximum  energies.  This  is  in  agreement  with  measurements.  Second, 
if  we  assume  that  the  emission  of  ions  lasts  for  200  nanoseconds,  not  an 
unreasonable  length  of  time,  we  find  that  the  corresponding  maximum 
energies  of  the  t  *’0  pulses  in  question  become  225  and  188  ev  respectively, 
in  good  agreement  with  the  measured  values  of  ion  energies  and  still  pre¬ 
serving  the  proper  difference  in  energy  between  the  two.  Third,  if  a 
finite  pulse  width  of  200  nanoseconds  is  assumed,  it  is  possible  that  a 
superposition  of  Maxwellians  for  that  time  interval  would  lead  to  measured 
distributions  such  as  these  of  Figures  II1-2  and  III- 3 .  Some  attempts  have 
been  made  to  simulate  such  finite  pulses  numerically,  as  reported  in  III  A  3; 
however,  more  work  needs  to  be  done  in  that  area  because  there  are  many 
adjustable  parameters. 

In  conclusion,  we  note  that  the  ion  energies  for  different  species  emitted 
from  a  given  target  are  approximately  the  same,  but  not  exactly  equal.  From 
results  obtained  so  far.,  it  appears  that  there  are  differences  between  the 
velocity  distributions  and  Maxwellian  distributions, 

2.  Velocity  Spectrum  of  Neutral  Molecule  Pulses 

The  analysis  of  the  velocity  distribution  can  also  be  applied  to  the  neutral 
particle  pulses  described  in  Section  I  to  determine  their  initial  energy 
distribution.  This  application  involves  more  uncertainties  than  in  t he  case 
of  the  ions  in  the  time  of -flight  instrument  since  the  particle  identities  are 
not  certain. 
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Let  us  carry  through  the  analysis  on  a  pulse  that  appears  typical.  The 

pulse  analyzed  is  shown  m  Figure  111  4.  This  data  was  obtained  wi’h  a 

2 

tungsten  target  and  a  laser  power  of  about  60  megawatts /cm  There  are 
two  different  particle  pulses,  This  icsult  is  ’vpieal  of  the  data  obtained 
with  Tungsten  tat  gets 


These  pulses  were  analyzed  according  to  the  mrthod  desci  ilx-tl  ear.u'i 
The  resulting  velocity  spec* ra  are  shown  m  Figures  111  5  ant)  ill  6.  Also- 
shown  in  each  figure  is  the  Maxwellian  distribution  for  the  same  average 
velocity,  and  normalized  to  the  same  number  of  particles  as  the 
experimental  Nnv..*.  The  results  of  Figure  11!  3  are  vcrv  similar  o  those 
descr  ibed  above  m  ‘he  analysis  of  the  ion  pulses.  The  veiucrc  dis*  nbucon 
appears  to  consist  o’’  a  du  cc’**d  componerr  and  a  tandem  component  As 
before,  the  Maxwe  i.an  is  not  a  good  fi»  The  observed  cut  vo  has  a  peak 
at  higher  value  ’ban  w,  and  fails  oft  mote  rapidly  on  *h<-  high  v<  .o<-  n  v  end 
Neat  ’he  wings,  as  ’he  distribution  upptoaches  /*.  t  o.  'he-  ohs<  i  v  •  <  1  <1;s'm 
but  ion  fa. is  bcT-w  ’he  Max  wo- Ilian.  The  appeal  anc»-  o’  tin  peak  at  a  \-io<  \ 

larger  ’ban  w  is  similar  *o  the  res  ul»  obtained  with  *  he  on  pu  isos  \V  do 
not  ex  pec*  ’he  distr  i  but  ion  To  go  to  zero  as  the  tens  did  since  7  he  eft  .<  *  v 
Y  is  zoic  a*  w  0. 

In  deriving  these  distributions,  \v<  have  ne  greeted  tie  fact  ’  hat  the  otfn  ler.cy 
y,  w:'h  which  a  high  voice  nv  atom  produces  elec* tens  at  the  cathode  o!  ’he 
muhipiier ,  is  a  func’ion  of  energy  L’ttle  da’u  on  this  efficiency  is  a\ui!ab'e 
inn  setne  estimates  ot  ns  pc s ible  behavior  ate  made  m  S**o» urn  I.  If  the 
elt  iciciu  y  vat  t e s  us  a  limcict)  of  atom  v«  Ion*  v  over  'in  range  of  m’erest,  n 
means  that  tin  eui  s  pt  •  s*  died  here  should  be  modified  bv  a  iae'oi  coil 
i  a  in  m  g  '  lie  e(  f  ie  iem  v  .  In  pa  r  t  ieu  !a  t  if  .  ttu  i  ■  uses  w  i  *  h  me  reus  mg  v**  !o<  ’  n  v 
cl  'lie  'iieulein  paiMeie,  'he  low  t  net  g  v  pot  '  ion  ol  the  d's’ t  dun  ion  should  b< 
me  t  *  usi  d  m  bt'  i  \  e  tc:  the  It?  gh  eno  f  g\  ;k  i  'u-n,  VV*  -  ir<  nd  ’  o  urn  s'  -,i't  -  las 
far*  oi  fur'll*  r  . 


Figure  III- 4  -  Typical  High  Speed  Molecule  Pulse  Obtained  by  Irradiation 
of  Tungsten  Target  with  Laser  Pulse  of  60  Mw/cra2.  Time  increases  from 
left  to  right. 
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Figure  III- G  -  Velocity  Distribution  N  (w)  Calculated  from  Second  Neutral 
Particle  Pulse  of  Figure  III-4.  Maxwellian  Fit  Included  for  Comparison. 
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Since  the  natures  of  the  particles  are  not  known  yet,  we  give  the  results  fnr 
the  energy  distributions  in  terms  of  the  atomic  number  A  of  the  particle. 

For  the  first  molecule  pulse  in  Figure  111-4,  we  find: 

Center  of  Mass  Energy  in  ev  =  39.0  x  A 

Random  Energy  -  one  component  -  in  ev  =  0.  91  x  A 

Random  Energy  -  three  components  -  in  ev  =  2.  7  x  A 

For  the  second  neutral  molecule  pulse  in  Figure  Ill  -  4  we  find: 

Center  of  Mass  Energy  in  ev  -  1.  57  x  A 

Random  Energy  three  components  -  in  ev  0.  45  x  A 

If  the  atomic  number  A  of  the  particle  were  known,  the  energy  would  be 
known;  eg.,  if  the  first  pulse  has  A  r  2  (molecular  hydrogen),  we  would 
obtain  78  ev  for  the  cenrer  of  mass  energy  and  5.  4  ev  for  the  random  enei  gy. 

We  should  note  that  this  is  one  sample  which  we  have  described  in  de  a-I 
here  Other  pulses  have  been  analyzed  and  yield  similar  resuK-  both.  fo>  •  re¬ 
shape  of  the  distributor.  ami  the  energy. 

These  results  win  be  used  below  in  Section  UTC  where  we  will  attempt  to 
deduce  the  nature  of  these  atom  pulses.  The  conclusions  of  this  subset-’ ion 
are  that  the  same  analysis  that  was  applied  to  the  ion  pulse  shapes  may  be 
applied  to  the  neutral  a*om  pulse  shapes.  The  results  indicate  a  la>  re¬ 
directed  energy  component  with  a  smaller  random  component.  Th<-  \  ■  o.  \ 
distributions  differ  somewhat  from  a  Maxwellian.  In  these  respects  'h- 
results  are  similar  to  those  obtained  in  the  analysis  of  the  ion  pulses.  Tin 
magnitudes  of  the  energy  components  cannot  be  deduced  here  because  of 
uncertainty  about  the  nature  of  the  particles,  but  only  the  ratio  r>f  energy  -o 
atomic  number. 


-43- 


3. 


Calculation  of  Pulse  Shapes  From  Inferred  Velocity  Distributions 


This  subsection  has  to  do  with  an  attempt  to  remove  the  6  function  assumption 
in  the  nalysis  described  in  Sections  111  A 1  and  IIIA2.  In  Appendix  1  it  is 
suggested  that  the  shape  of  the  pulse  J  <t)  be  calculated  using  the  assumptions 
of  a  square  emission  pulse  in  time  and  a  Maxwellian  distribution  for  the 
velocity.  This  procedure  involves  evaluating  the  integral  in  Equation  .23'  of 
Appendix  1.  As  is  pointed  out  there,  even  with  these  fairly  restrictive 
assumptions,  the  integral  cannot  be  evaluated  analytically. 


Accordingly  a  computer  program  was  written  to  evaluate  this  integral 
numerically.  The  case  chosen  was  that  of  mass  39  ions  with  a  1000  volt 
accelerating  potential.  The  shape  of  the  pulse  was  calculated  by  evaluating 
the  integral  using  Equation  (23)  described  m  Appendix  1.  The  function  r  , 
defined  in  Appendix  1  as  the  density  of  ions  per  unit  initial  speed  w  in  the  x 
direction,  is  given  for  a  Maxwellian  by: 


c  io.  w.  t)  =  F't)  exp 


m(w-w) 

2kT(t) 


(111-  7) 


Here  w  is  the  average  velocity  of  the  distribution,  corresponding  to  the 
directed  energy  compcnent;  it  is  assigned  as  an  independent  variable  so  that 
the  peak  of  the  pulse  occurs  at  the  same  time  as  that  experimentally 
observed.  The  following  assumptions  are  made:  (1)  The  random  component 
is  a  Maxwellian  characterized  by  a  temperature.  (2/  the  temperature  is  a 
given  function  of  time;  (3)  The  function  F(?)  gives  the  shape  of  the  emission 
pulse  and  is  also  known.  Table  III- 2  lists  the  various  cases  which  were 
calculated.  Equation  <23 *  strictly  applies  only  to  those  cases  in  which  the 
emission  pulse  is  flat  and  the  temperature  is  constant.  Expressions  some¬ 
what  different  from  the  integral  in  Equation  (23)  must  be  evaluated  for  the 
cases  in  the  table  which  do  not  comply  with  these  restrictions. 
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SETS  OF  EMISSION  PULSE  PARAMETERS  USED  IN  THE  INTEGRATION  OF 

EQUATION  (23)  OF  APPENDIX  I. 


CASE 

W 

(10hC  M/SEC) 

F  (0 

T  (0 

DURA  HUN  OF 

PUL  St  (NS) 

SHAPE  or 

PULSL 

SHAPE 

PEAK  TEMPERATURE 

"K 

1 

r.5 

200 

FLAT 

f  LAT 

105 

2 

2.5 

200 

FLAT 

FLAT 

2X105 

3 

4.0 

200 

FLAT 

FLAT 

106 

4 

4.0 

20C 

FLAT 

FLAT 

105 

5 

3.0 

200 

FLAT 

FLAT 

105 

6 

3.0 

200 

TRIANGULAR 
(PEAK  AT  0) 

FLAT 

105 

7 

3.0 

200 

TRIANGULAR 
(PEAK  AT  0) 

FLAT 

2X105 

8 

3.0 

200 

FLAT 

FLAT 

2X1 r  3 

9 

2.5 

100 

FLAT 

FLAT 

105 

10 

2.5 

200 

FLAT 

FLAT 

3X103 

11 

2.5 

200 

FLAT 

FLAT 

6X103 

12 

2.5 

200 

FLAT 

FLAT 

12X103 

13 

3.5 

200 

FLAT 

FLAT 

104 

14 

3.5 

200 

FLAT 

FLAT 

103 

15 

3.5 

200 

TRIANGULAR 
(PEAK  AT  0) 

FLAT 

103 

16 

3.5 

100 

FLAT 

FLAT 

103 

J'/ 

3.5 

50 

FLAT 

FLAT 

105 

18 

5.5 

100 

FLAT 

FLAT 

2X105 

19 

3.5 

200 

FLAT 

FLAT 

2X103 

20 

3.5 

200 

FLAT 

TRIANGULAR 
(PEAK  AT  0) 

103 

21 

3.5 

1000 

FLAT 

FLAT 

103 

22 

3.5 

1000 

TRIANGULAR 
(PEAK  AT  0) 

FLAT 

103 

23 

3.5 

1000 

TRIANGULAR 
(PEAK  AT  0) 

FLAT 

5X104 

24 

3.5 

1000 

TRIANGULAR 
(PEAK  AT  0) 

TRIANGULAR 
(PEAK  AT  0) 

103 

25 

3.5 

1000 

FLAT 

TRIANGULAR 
(PEAK  AT  0) 

103 
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This  computation  essentially  does  in  reverse  the  analysis  done  in  Sections  T11A 
and  LTA3.  There  we  started  with  an  observed  pulse  shape  and  obtained  a 
velocity  distribution.  Now  we  start  with  an  assumed  velocity  distribution 
and  calculate  a  pulse  shape. 

6 

We  find  that  a  choice  of  w  '  3.  5  x  10  cm/sec  gives  the  best  fit  for  the  peak 
of  the  calculated  data  as  compared  to  the  experimentally  observed  pulses. 
Most  of  the  cases  in  Table  III- 2  were  run  with  this  value  of  w  . 

From  the  computer  runs,  we  found  that  the  calculated  pulse  shapes  could  not 
be  varied  appreciably  by  any  of  the  choices  of  parameters  listed  in  Table 
III- 2 .  As  was  pointed  out  before,  care  must  be  exercised  in  choosing  the 
mesh  size  for  the  numerical  integration.  The  first  few  runs  on  the  computer 
were  made  using  too  large  a  mesh  which  gave  spurious  results  indicating 
considerable  changes  of  the  pulse  shape  with  the  variable  parameters.  Later 
work  with  a  smaller  mesh  is  described  in  this  reoor:. 

Some  of  the  results  of  the  computer  runs  are  presented  in  Figures  111  7,  8, 

9  and  10.  In  these  figures  the  calculated  pulse  shapes  are  shown  for  sevei  al 

6 

of  the  cases  listed  in  Table  III- 2,  all  of  them  with  w  -  3.5  x  10  cm/sec. 

The  vortical  scales  are  in  arbitrary  units  they  are  not  even  the  same  scale 
for  all  cases.  The  important  consideration  is  the  relative  shape  of  the 
pulses;  e.  g,  ,  the  time  of  the  peak,  the  width,  the  relative  slope  of  the 
rising  and  decaying  portions. 

In  Figure  111-7  we  see  the  effect  of  changing  the  duration  of  the  emission  Fit) 
from  50  to  1000  nanosecond,  with  Fit)  assumed  flat  for  all  the  cases.  As 
the  length  of  F't)  increases,  the  width  of  J(t)  increases,  as  we  would  expect, 
and  the  peak  value  of  J(t)  comes  later,  as  we  also  expect.  However,  the 
shape  of  the  pulse  Jit)  does  not  change  dramatically  it  remains  roughly  a 
symmetrical  triangular  pulse.  If  the  velocity  distribution  obtained  in 


Figure  III- 7  -  Detector  Pulse  Shapes  Claculated  Using  a  Constant  Emission  F(t)  and 
Constant  Temperature.  The  Different  Curves  Represent  Different  Widths  of  F(t), 
from  50  to  1000  Nanoseconds.  The  Number  on  Each  Identifies  the  Set  of  Parameters 
in  Table  III- 2  Corresponding  to  that  Curve. 


-47- 


(ARBITRARY  UNITS) 


Figure  III -8  -  Detector  Pulse  Shapes  Calculated  Using  a  Constant  Emission  F(t)  and 

&  4  5 

Constant  Temperature.  Different  Values  of  Temperature,  from  10  to  2  x  10"  °K. 

The  Number  on  Each  Identifies  a  Set  of  Parameters  in  Table  III -2  Corresponding  to 
that  Curve. 
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Figure  III  -  9  -  Detector  Pulse  Shape  Calculated  Using  a  Constant  Emission  F(t)  and 
Time-dependent  Temperature  T(t).  The  Two  Curves  Represent  Different  Shapes 
for  T(t).  The  Number  on  each  Identifies  a  Set  of  Parameters  in  Table  III  - 2  Corres 
ponding  to  that  Curve. 
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Figure  III- 10  -  Detector  Pulse  Shape  Calculated  Using  a  Time-Varying  Emission  F(t) 
and  Constant  Temperature.  The  Two  Curves  Represent  Different  Shapes  for  F(t). 

The  Number  on  Each  Identifies  a  Set  of  Parameters  in  Table  III-2  Corresponding  to 
that  Curve. 


Section  III  A  2  using  a  6  function  approximation  were  Maxwellian,  we  would 
expect,  in  these  calculations,  that  as  we  shorten  the  length  of  the  emission 
F(t),  we  should  approach  a  pulse  shape  similar  to  that  observed  in  the 
experiments.  But  this  is  not  the  case.  We  shall  return  to  this  point  later. 

In  Figure  III- 8  we  see  the  effect  of  increasing  the  temperature  T(t) 

4  5 

characterizing  the  Maxwellian  from  10  to  2  x  10  degrees  Kelvin,  for  flat 
F(t).  At  the  lower  end  of  this  range.,  we  find  pulses  too  short  to  agree  with 
experiment,  and  at  the  upper  end  pulses  too  long,  but  again  the  pulse  shape 
does  not  change  significantly. 

In  Figure  III- 9,  we  see  the  effect  of  changing  the  shape  of  T(t*,  without 
changing  the  peak  value,  again  there  is  no  great  effect  on  the  shape  of  J<t>. 

In  Figure  III- 10,  we  see  the  effect  of  changing  the  shape  of  F(t)  here  a 
detailed  comparison  of  the  curves  shows  that  a  triangular  F C t )  gives  a  J<t> 
that  rises  slightly  more  slowly  than  that  given  by  a  flat  F<tl,  but  the  effect 
is  not  great. 

As  an  extension  of  the  ideas  discussed  in  connection  with  Figure  III- 7  and 
as  a  further  confirmation  of  the  conclusions  reached  above  we  have  taken  F(t) 
to  be  a  delta  function  and  carried  out  the  integration.  The  results  are 
shown  in  Figure  III- 11.  Curve  1  was  obtained  using  the  velocity  distribution 
derived  in  the  first  use  of  the  analysis  described  in  Appendix  I,  which  was 
similar  to  a  Maxwellian;  as  noted  before  this  was  obtained  with  too  coarse  a 
mesh.  Curve  2  was  obtained  using  a  Maxwellian  velocity  distribution  with  a 
temperature  of  70,  000  degrees  Kelvin,  and  curve  3  was  obtained  using  an 
asymmetric  velocity  distribution  with  a  peak  near  the  high  velocity  end  of  the 
distribution,  similar  to  the  types  found  in  Section  III  A  1.  Curves  1  and  2 
are  similar  to  those  obtained  above.  The  results  of  curve  l.in  particular, show 
that  the  earliest  results,  which  cannot  indicate  a  velocity  distribution  similar 
to  a  Maxwellian,  are  not  exactly  right,  since  running  the  analysis  backwards 
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Figure  III- 11  -  Detector  Pulse  Shapes  Calculated  Using  a  Delta  F'unction  for  F (t ) 
and  Constant  Temperature.  Curve  1  is  from  a  Quasi-Maxwellian,  as  described 
in  text;  curve  2  is  from  a  Maxwellian  at  70,  000°K;  curve  3  is  from  an 
asymetric  velocity  distribution  peaked  near  the  hig'  /eloeity  end. 
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with  the  distribution  obtained  does  not  yield  the  original  pulse  shape  J(t). 
Curve  2  indicates  that  shortening  the  length  of  the  emission  F(t)  toward 
zero  does  not  modify  the  pulse  shape,  as  we  said  before  in  connection  with 
Figure  III- 7.  Curve  3  indicates  that  the  experimentally  observed  pulse  shape 
may  be  obtained  with  velocity  distributions  of  the  type  found  in  Section  III  A  1 
slightly  asymmetric  with  peaks  near  the  high  velocity  end. 

Thus,  this  analysis  so  far  indicates  that  the  calculated  pulse  shape  is  not 
particularly  sensitive  to  variation  of  the  available  parameters,  so  long  as 
a  Maxwellian  velocity  distribution  is  assumed.  An  exception  to  this  state¬ 
ment  may  occur  when  F(t)  varies  over  a  long  emission  time  and  T(t)  is  low. 
Further  work  is  in  progress  on  this  point.  The  preliminary  results  support 
the  results  of  Section  111  A  1  where  we  found  distributions  that  depart  from 
Maxwellian  form. 

B.  EXTENSION  OF  CALCULATIONS  ON  EFFECT  OF  INVERSE 
BREMSSTRAHLUNG 

In  the  last  semiannual  report*  we  described  calculations  on  the  possible  effect 
of  heating  of  the  blow-off  material  by  free-free  transitions  of  the  electrons 
in  the  path  of  the  laser  radiation,  and  came  to  the  conclusion  that  electron 
densities  were  too  low  to  explain  the  high  (200  ev)  ion  energies  observed. 

In  view  of  the  results  described  above,  in  which  it  appears  that  the  major 
portion  of  this  energy  is  a  directed  energy,  it  appears  worthwhile  to  re¬ 
examine  this  result  to  see  if  the  random  component,  an  order  of  magnitude 
lower,  could  possibly  be  explained  by  this  mechanism.  This  would  not  serve, 
to  explain  the  directed  component,  but  a  random  component  of  the  order  of 
10-20  ev  does  appear  to  be  similar  to  a  thermal  energy.  This  random  energy 
also  appears  to  be  compatible  with  the  temperatures  deduced  for  the  blow-off 
material  using  methods  cf  optical  spectroscopy. 
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We  consider  two  possible  means  by  which  the  effects  considered  in  the  semi¬ 
annual  report  could  be  enhanced;  i.e.,  increased  absorption  of  light  because 
of  sodium  and  potassium  ions  emitted  by  the  surface,  and  the  possibility  of 
inverse  Bremsstrahlung  with  the  electron  scattering  in  an  interaction  with  a 
neutral  particle  rather  than  an  ion.  The  first  involves  the  presence  of  alkali 
metal  ions  in  the  blow-off  material;  these  materials  which  would  be  present 
as  ions  and  not  as  neutral  atoms  would  increase  the  ion  density  over  the 
thermal  ion  density  of  the  dominantly  un-ionized  blow-off  material, and  thus 
increase  the  absorption  cross  section  for  the  light  in  the  critical  early 
stages  where  the  heating  is  just  beginning,  and  generally  is  not  proceeding 
rapidly  because  of  low  charged  particle  densities. 

This  consideration  involves  a  feedback  from  the  empirical  results  into  the 
theory.  We  would  be  using  our  knowledge  of  what  the  blow-off  material 
actually  consists  of.  The  theoretical  model  then  deals  with  a  case  more 
like  what  is  actually  realized. 

Accdrdingly  the  computational  model  described  in  the  semi  annual  report  has 

been  modified  by  adding  to  the  ion  density  a  term  corresponding  to  an  ion 

density  of  alkali  metal  ions.  This  added  ion  density  is  not  increased  as  the 

temperature  increases  because  we  consider  these  materials  to  be  essentially 

fully  ionized  when  they  are  emitted  from  the  surface.  We  estimate,  from  the 

12  3 

experimental  results,  an  alkali  metal  ion  density  of  about  10  /cm  in  the 
early  stages  of  the  development  of  the  particle  emission. 

The  results  of  these  calculations  appear  in  Figures  111-12,  13  and  14. 

i 

Figure  III- 12  is  identical  with  Figure  IV- 5  of  the  semi-annual  report*  and 
represents  the  case  without  the  added  alkali  ions  for  a  low  initial  neutral 
particle  density.  Figure  III- 13  shows  the  effect  of  adding  the  additional 
ions;  the  temperature  increases  much  more  rapidly  for  a  constant  initial 
electron  density,  and  the  critical  electron  density  for  runaway  heating 


-54- 


MMm,  ,uy *"■— 
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Figure  III  - 1 2  -  Calculated  Plasma  Heating  by  Inverse  Bremsstrahlung. 
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RATIO  OF  ION  DENSITY  WITH  ABSORPTION  TO  ION  DENSITY  WITHOUT  ABSORPTION 


Figure  III  - 1 3  -  Calculated  Plasma  Heating  by  Inverse  Bremsstrahlung  in  the  Presence 
of  Fully  Ionized  Sodium 
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RATIO  OF  ION  OENSITY  WITH  ABSORPTION  TO  ION  DENSITY  WITHOUT  ABSORPTION 


Figure  III  - 1 4  -  Effect  of  Increased  Neutral  Particle  Density  on  Plasma  Heating  by 
Inverse  Bremsstrahlung 
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RATIO  OF  ION  DENSITY  WITH  ABSORPTION  TO  (ON  DENSITY  WITHOUT  ABSORPTION 


during  the  period  of  the  laser  pulse  is  reduced  by  a  factor  of  five.  Figure 
III- 14  shows  the  results  at  a  higher  initial  neutral  particle  density.  The 
curves  are  essentially  the  same  for  10  alkali  ions /cm  as  for  0.  In  this 
case,  the  addition  of  the  alkali  ions  makes  no  difference  since  the  thermal 
ionization  yields  an  ion  density  greater  than  that  of  the  ions.  Comparison 
of  Figures  III  - 12  and  III-  13  shows  that  the  critical  electron  density  is  reduced 
at  higher  initial  neutral  particle  densities  although  the  variation  is  rather  slow. 
The  conclusion  is  that  the  presence  of  the  alkali  metal  ions  should  indeed  make 
the  heating  proceed  more  rapidly,  especially  for  relatively  low  particle 
densities,  and  that  the  restrictions  on  electron  density  i'or  heating  are  eased 
somewhat.  These  changes  by  themselves,  however,  are  not  sufficient  to 
reverse  the  conclusions  in  the  semi-annual  report. 

Now  let  us  consider  the  possibility  of  inverse  Bremsstrahlung  proceeding  by 
interaction  of  the  electron  with  a  neutral  particle  rather  than  with  an  ion. 

7 

There  have  been  seme  recent  results  on  the  cross  section  ior  this  process. 

The  cross  section  i-s  smaller  than  for  the  process  involving  an  ion,  but  this 
may  be  compensated  by  the  much  larger  density  of  neutral  particles  which  is 
available  in  the  early  stages  of  the  process.  When  we  also  consider  that  we 
may  be  trying  to  explain  a  thermal  energy  of  only  20ev  rather  than  200  ev  as  we 
were  before,  this  additional  process  may  be  helpful.  The  analysis  is  now  in  the 
process  of  being  carried  out,  but  no  results  are  available  as  yet. 

The  absorption  coefficient  k^  for  light  in  the  slightly  ionized  blow-off  material 
by  the  process  in  which  the  electron  makes  a  free -free  transition  in  the  field  of 
an  ion  was  given  in  the  semi-annual  report; 


k  (ion) 
v 
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for  the  ruby  laser  frequency.  Here  n  and  n.  are  electron  and  ion  densities, 

T  is  the  temperature  of  the  material,  v  is  the  frequency  of  the  light,  and  in  A 
accounts  for  collective  effects  and  is  approximately  equal  to  10.  The  symbols 

? .  .  '  j 

m,  e,  c,  and  k  have  their  conventional  meanings  of  electron  mass,  electron 
charge,  velocity  of  light,  and  Boltzmanms  constant. 


For  the  free-free  transition  of  an  electron  near  a  neutral  hydrogen  atom,  the 

7 

absorption  coefficient  has  recently  been  evaluated  as 


(neutral) 


4.&V'2 

3  3/2 

m  c  v 


n  n  =  7.  79  x 
n  e 


n 


n 


for  the  ruby  laser  frequency. 

Here  nn  is  the  density  of  neutral  hydrogen  atoms,  and  ti  is  Planck's  constant 
divided  by  2tt  .  This  result  is  strictly  valid  only  for  atomic  hydrogen; 

8 

corresponding  results  for  other  gases  are  noT  available.  One  estimate 

for  molecular  hydrogen  indicates  an  absorption  coefficient  smaller  by  a 

9 

factor  of  between  two  and  four.  An  estimate  for  argon  indicates  a  factor 
smaller  by  a  factor  of  four.  However,  if  we  assume  that  the  absorption  cross 
section  for  the  other  materials  emitted  as  neutral  gases  (  such  as  CO  and  COg 
which  are  known  to  be  produced  copiously  under  our  experimental  conditions)  is 
comparable  to  that  of  hydrogen,  we  are  led  to  a  large  value  of  the  absorption 
coefficient  by  this  process;  indeed,  to  a  larger  value  than  that  due  to  the 
absorption  in  the  field  of  an  ion.  The  ratio  of  the  absorption  coefficients  due 
to  the  two  processes  would  be: 


k  (neutral) 
v  _ 

k  (ion) 
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Typically,  for  emission  of  the  usual  gases  from  a  surface  at  a  few  thousand 
degrees  Kelvin,  we  would  expect 

n./n  10  ®  to  10 
in 

so  that  the  ratio  defined  above  could  be  much  greater  than  one.  Thus, 
absorption  in  collisions  of  electrons  with  neutral  atoms  and  molecules  could 
easily  dominate  over  absorption  in  collisions  with  ions.  If  this  were  to  occur, 
the  critical  electron  density  at  which  heating  occurs  rapidly  could  be  lower 
than  calculated  in  the  semi-annual  report. 

The  computations  based  on  these  considerations  are  similar  to  those  based  on 
the  absorption  in  collisions  with  ions  described  in  the  semi-annual  report. 

The  experimental  studies  of  neutral  gas  desorption  allow  realistic  estimates 
of  neutral  particle  density  as  a  function  of  time  near  the  surface.  The  result 
of  carrying  out  this  procedure  with  the  larger  absorption  coefficients  may 
predict  sufficient  heating  to  explain  the  random  component  of  the  velocity 
distributions  described  earlier,  but  they  are  unlikely  to  explain  the  directed 
component. 

C.  TENTATIVE  IDENTIFICATION  OF  HIGH  SPEED  NEUTRAL  MOLECULES 
EMITTED  FROM  TUNGSTEN 

Let  us  now  consider  the  nature  of  the  particles  contained  in  the  neutral  pulses. 
We  use  the  results  obtained  in  Section  III  A  2  on  the  data  of  Figure  III- 4. 

In  that  section  we  obtained  results  for  the  energy  in  terms  of  the  mass  number 
A  for  both  the  directed  and  random  components. 

If  we  were  to  assume  that  the  first  neutral  particle  pulse  was  CO,  an 
assumption  that  might  seem  plausible  since  CO  is  the  most  abundant  com¬ 
ponent  observed  in  the  gas  phase  emitted  from  tungsten,  we  note  that  the 
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directed  energy  would  be  approximately  12  70  ev.  This  does  not  seem 
reasonable  since  the  energy  would  be  much  higher  than  that  found  in  the 
measurements  on  ions  in  the  time  of  flight  instrument.  Moreover,  let  us 
consider  the  ratio  of  the  masses  represented  by  the  two  pulses.  If  we 
assume  that  the  two  peaks  represent  approximately  the  same  energy,  we 
find 


1 

2  mi 


where  t.  is  -he  time  of  arrival  of  mass  m.  (measured  at  the  peak  of  the  pulse) 
and  d  is  the  distance  from  target  to  detector.  Then 


mi  \h  I 

At  the  peaks  of  the  pulse  tj  “  6.04us,  t2  ~  22.  2us,  leading  to  mg/m^  13.5. 

If  the  first  pulse  were  CO,  the  second  would  then  have  a  mass  number  of  almost 
380,  an  unreasonably  high  value.  If  we  were  to  assume  the  second  pulse  to  be 
CO,  the  first  would  then  be  mass  2  which  would  correspond  to  H2>  This  appears 
to  be  a  reasonable  component  to  be  emitted.  If  the  identities  of. these  two 
pulses  are  H2  and  CO,  the  energies  are,  according  tc  the  results  of  Section 
III  A  2,  78  and  44  ev  for  the  directed  components,  and  5.  4  and  12.  6  ev  for 
the  random  components. 


These  energies  are  not  in  perfect  agreement,  they  violate  the  assumption  of 
equal  energy  for  the  ’wo  pulses  which  we  advanced  above.  They  are  also 
somewha*  lower  than  the  energies  found  in  the  ton  emission  work,  but  the  78 
ev  for  H9,  at  any  rate,  differs  only  by  a  factor  of  2  from  the  ion  energies. 
These  energies,  at  any  rate,  appear  tc  be  more  plaut  ble  than  the  1270  ev 
found  by  postulating  that  the  first  pulse  is  CO  Ther*  are  no  other  easonable 
mass  assignments  for  »hese  two  peaks  d “'•  r  ing  by  a  factor  of  14.  In 
midi*  on,  as  we  noted  in  Section  1,  there  is  evidence  for  identifying  the  first 
neutral  molecule  peak  frem  titanium  n t  hydrogen.  Clearly  a  considerable 
amount  of  additional  work  is  required  to  explain  these  results. 


SUMMARY 


The  work  described  in  this  report  deals  with  the  measurement  of  particle 
emission  from  metal  surfaces  irradiated  with  light  flux  densities  of  50-70 
Mw/cm  ,  as  well  as  calculations  of  the  velocity  distribution  of  the  particles 
from  the  measured  pulse  shapes.  Measurements  in  the  quadrupole  mass 
spectrometer  have  shown  that,  besides  the  thermal  gas  desorption,  the  neutral 
particle  emission  contains  a  photon  pulse  and  several  high  velocity  non-zero- 
mass  pulses.  Tentative  assignments  of  masses  and  energies  of  the  high  speed 
neutral  pulses  are  made. 

The  electric  field  dependence  of  the  time  of  flight  of  ions  emitted  from 
tungsten  and  platinum  surfaces  were  measured, and  it  was  concluded  that 
they  behave  as  independent  particles  in  their  response  to  electric  fields.  The 
ion  mass  spectra  obtained  from  platinum  were  very  similar  to  those  from 
tungsten  with  respect  to  both  mass  and  energies  observed.  Other  measure¬ 
ments  on  platinum  suggest  that  an  interaction  of  the  blow-off  particles  with 
the  laser  beam  is  required  for  the  production  and  acceleration  of  ions. 

Analysis  of  *he  measured  pulse  shapes  for  both  ions  and  neutrals  yields 
velocity  distributions  with  random  energies  of  a  few  tens  of  electron  volts 
superimposed  on  a  directed  energy  of  hundreds  of  electron  volts.  Cal¬ 
culations  of  heating  of  the  blow-off  material  by  inverse  Bremsstrahlung 
have  been  extended  to  the  case  of  a  multi -component  system,  with  one  of 
the  components  being  a  fully  ionized  alkali  metal. 
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APPENDIX  I 


ION  VELOCITY  DISTRIBUTION  FROM  MEASURED  MASS  PULSES 


The  material  m  this  appendix  consists  of  unpublished  work  by  Frank  J.  'lien 
of  Ballistic  Research  Laboratories.  It  is  included  here  because  we  draw 
heavily  on  its  contents,  and  the  material  is  not  available  elsewhere.  We  have 
edited  and  condensed  the  original,  retaining  only  those  parts  required  for  use  in 
this  report. 


Suppose  the  particles  vaporized  by  the  laser  pulse,  whether  neutral  or  ionized 
initially,  form  a  dense  plasma  close  to  the  target  surface.  Many  interactions 
take  place  here;  the  particles  also  undergo  collisions  with  the  target  surface. 
The  result  of  all  this  is  assumed  to  be:  ions,  electrons,  and  neutral  atoms 
(and  in  some  cases  more  complex  species)  attain  random  velocities  superposed 
upon  a  dir  ected  component  of  velocity  normal  to  the  surface.  We  might  expect 
for  each  species  a  Maxwellian  distribution  about  the  directed  component,  but 
this  must  be  decided  by  the  experimental  results. 

The  ions  after  leaving  the  placm.  i.e,  after  getting  far  enough  from  the 
target  surface  so  that  the  density  is  low  and  collisions  in  which  significant 
momentum  is  transferred  are  infrequent,  are  accelerated  in  an  electric  field 
and  ultimately  those  traveling  in  a  very  small  solid  angle  about  the  normal 
to  the  target  surface  are  detected.  Under  these  conditions  we  might  expect 
that  the  mere  energy  the  ions  have,  the  higher  will  be  both  their  random 
velocities  and  the  superposed  directed  velocity  normal  to  the  target  surface. 

Before  considering  a  distribution  of  ion  velocities,  consider  the  detection  of  a 
single  pulse  of  ions  with  velocity  w  normal  to  the  target  surface.  The  ions 
travel  one  dimensionally  in  a  region  consisting  of  an  arrangement  of  grids  with 
potentials  which  are  completely  arbitrary  except  for  the  restriction  that  the 
ions  not  be  reversed  in  direction,  i.e.,  the  ions  are  allowed  to  reach  the 
detector  at  the  end  of  the  region.  Let  p,  v,  J  be  the  ion  density,  velocity, 
and  current  and  V  the  field  potential.  Then  the  governing  equations  are: 

d  J  _  _  do  /  \ 

ax  at  *  ' 

J  “  p  v  (2 ) 

-  m  v2  (x)  2  =  e  V(x)  -  V(0)  (3) 

2  1  L 

where  e  is  the  Ion  charge. 
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The  "initial  current"  is  J(0,t)  =  p(0,t)  w.  Now  the  ions  emitted  at  time  t 
(i.e.,  leave  the  plasma  without  further 'collisions)  reach  a  point  x  at  a  time 
t  +  t  where 


(4) 


where  v(x)  =  ~\J  —  £  V(x)  -  V(0)  j  +  w 


Nov/  the  ions  emitted  in  time  dt  at  x  =  0  fill  a  tube  of  length  w  dt  while  at 

time  t(x)  later  these  same  ions  fill  a  space  v(x)  dt.  But  the  number 

crossing  a  plane  at  (x,  t)  must  be  the  same  as  the  number  emitted  at  {0,  t-  r ). 

J  (x,  t)  =  p  (x,  t)  v  (x)  =  o(o,  1-t)  w  -  J  (o,  t-~)  (6) 


It  is  readily  verified  that  this  prescription  satisfies  the  continuity  equation: 

SJ  _  3  o  (o,  t-r)  3(t-^) 

* -  -  w  -  x  - 

3 x  o  (t-T)  ex 


-  w 


But  T  (x) 


3o(o,  t-T  ) 

3(t-,r) 

dx 
v  (x) 


-  St 
d  x 


Thus  -^L 
Sx 


v(x) 


j  W  x  Sp  (0,  t-T) 
^v(x)  O(t-T) 
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Now  from  equation  (5): 

-  ap  =  j 
at  \ 

=  !  ~w  \  dp(0,  t-T ) 

^  v(x)  j  a  (t  -  t  ) 

So 

ax  a  t 

Thus  the  governing  equations  are  satisfied  no  matter  how  many  grids  there  are, 
where  they  are,  or  what  the  potentials  are,  provided  only  the  potentials  are  such 
as  to  allow  the  ions  to  reach  the  detector.  This  condition  implies  v(x)  >  0  and 
t  (x)  >  0. 

To  recapitulate,  the  solution  is 

J(x,  t)  =  p(x,  t)  v(x)  =  o  (0,  t-T  )  w  =  J  (0,  t-T  ) 

dx 
v(x) 


v(x)  = 

Since  it  is  assumed  that: 

p  (0,  t)  =  0  for  t  <  0, 

J  (x,  0)  =  0. 

This  solution  has  the  important  property  that  the  current  pulse  J(x,  t)  suffers 
no  distortion  at  any  point  in  the  flow  field;  the  (retardation)  time  t  (x)  depends 
upon  the  grid  and  potential  arrangement,  but  the  current  pulse  J(x,  t)  =  J (0,  t-T  ). 


-  W  ?P  (0,  t-T  )  a  (t-T  ) 

-  -  x  - 

v(x)/  a  (t-r)  at 
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Consider  now  a  case  in  which  the  ions  leave  the  plasma  with  some  distribution 
of  speeds  in  the  x  direction  with  the  same  other  conditions  as  previously. 

(The  ions  with  significant  y  and  z  components  will  not  be  detected). 

Let  w  8  speed  in  x  direction  ion  has  upon  escaping  from  plasma. 

Let  p  (0,  Wjt)  8  density  per  unit  initial  speed. 

J(0,  w, t)  =  current  per  unit  initial  speed. 

Let  p  (x,  v,  t)  =  density  per  unit  ion  speed  v. 

J  (x,  v,t)  =  current  per  unit  ion  speed  v. 


oo  ® 

ThenJ(0,t)  =  /  p(0,w,t)wdw  *  I 

-4=0  *4 


J(0,  w,  t)  dw 


00  00 

J  (x,  t)  =  J  p  (x,  v,  t)  v  (x)  dv  =  J  J  (x,  V,  t)  dv 


0 


0 


(8) 

(9) 


By  analogy  with  the  solution  for  the  constant  velocity  case  it  then  appears  that 
the  solution  for  this  case  should  be  given  by  the  following  prescription: 


00 


J(x,  t)  =  / 
JQ 

D  [0, 

i  L 

w,  t,  -  T  (x,  w)  "J  w  d  w 

X 

T  (x,  w)  = 

f 

dx 

Jo 

V  (x.  w) 

V  (x,  w)  = 

r  V(x)  -  V(0)  ]  +  w2 

V  m 

L  J 

Again  p  (0,  w,  t) 

*  0  for  t 

^  o, 

and  v  (x,  w)  >  0  , 

T  (x,  w)  >  0, 

for  all  detected  ions.  Then  J  (x#  0)  ■  0  as  before. 
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Again  it  may  be  shown  that  this  physically  obtained  solution  satisfies  the 
continuity  equation  by  differentiating  J  (x,  t)  as  required.  It  is  helpful  to  note 
that  vdv  =  wdw  in  carrying  out  the  procedure.  The  procedure  is  the  same  as 
previously  and  will  not  be  repeated. 

A.  DELTA  FUNCTION  APPROXIMATION  FOR  EMISSION  PULSE 

The  solution  given  is  mathematically  correct  for  any  arbitrarily  prescribed 
function  J  (0,  w,  t).  However,  the  physical  problem  is  to  find  J  (0, w,  t)  from 
the  experimentally  determined  J(x,  t).  This  cannot  be  done  uniquely;  some 
assumption  has  to  be  made  about  J(0,  w,  t).  Since  the  laser  pulse  is  short  and 
the  laser-induced  plasma  has  a  fleeting  existence,  the  most  obvious  assumption 
is  that  J(0,  w,  t)  is  a  delta  function  in  time.  This  appears  reasonable  when  the 
time  spread  in  the  mass  peak  is  long  compared  with  the  laser  pulse  time. 

(The  lighter  mass  peaks  could  be  spread  out  somewhat  in  time  by  lowering 
the  accelerating  potential;  this  procedure  is  limited  by  the  rather  high  constant 
speed  of  the  plasma  normal  to  the  target  surface  on  which  the  random  motions 
are  superimposed,  and  by  the  decrease  of  solid  angle  in  which  emitted  ions  are 
detected). 

B’or  cases  in  which  a  delta  function  in  time  cannot  be  assumed,  one  no  longer 
has  a  one-to-one  correspondence  between  escape  velocity  w  and  time  of 
detection  t.  This  correspondence  leads  to  a  considerable  simplification  in  the 
nature  of  the  solution  given  by  Eq.  (10).  We  consider  the  case  first  therefore 
in  which  the  ions  escape  the  plasma  as  a  6  function  in  time. 

For  this  case  it  is  simpler,  and  physically  more  satisfying,  to  proceed  directly 
rather  than  to  use  the  solution  for  the  more  general  case.  First  note  that 
almost  all  of  the  time  the  ions  spend  in  getting  from  the  plasma  to  the  detector 
is  spent  in  the  drift  tube  where  each  ion  has  a  constant  velocity.  This  constant 
velocity  is  determined  by 

-  m  w2  +  e  [v  (xd)  -  V  (0)  ]  =  -  m  v2  (11) 

2  2 

where  v  is  the  velocity  in  the  drift  tube  and  V(x,-|)  =  potential  throughout  drift 
tube,  w  is  the  ion  velocity  upon  escaping  the  plasma  and  is  distributed  in  some 
way  about  some  mean  plasma  speed  away  from  the  target  surface.  Since  the 
detector  position  is  fixed,  we  replace  J  (x,  t)  by  J(t),  the  detected  current. 

Let  N(w)  =  the  number  of  ions  escaping  plasma  per  unit  escape  velocity  per 

square  centimeter  during  the  dbltai  function  escape  time. 

00 

Let  N  =  /n(w)  dw  (12) 
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In  the  previous  notation 
N  =  /  J  (0,  w,  t)  dw  6(t)  dt  = 


=  Jj  (0,  w( 


J  J  (0,  w. 


0)  dw 


Then 


N(w)  dw  =  J  (t)  dt 


(13) 

(14) 


where  t  and  w  are  refuted  since  t  is  the  time  spent  in  the  drift  tube  by  an  ion 
of  initial  escape  speed  w.  In  this  section  we  make  the  approximation  that  the 
time  not  spent  in  the  drift  tube  doesn't  contribute  to  the  spreading  of  the  mass 
peak;  if  the  time  of  detection  is  measured  with  respect  to  the  laser  pulse  zero 
time,  the  time  not  spent  in  the  drift  tube  is  a  constant  to  be  subtracted  from  the 
total  time  in  using  the  relationship  between  t  and  w. 

We  thus  have 

t  =  length  of  drift  tube  Xd 

velocity  in  drift  tube  v 


or 


and 


dt 

dw 


xdw 


3 

v 


xdw 


Then 


N(w) 


J(t) 


w 

xd 

l\x  d 

with  w  =  -v  /[ — 

2 

-  2e 

V(xd)  -  V(0) 

V  \  t 

m 

m  * 

(16) 


(17) 


(18) 


(19) 


Equations  (18)  and  (19)  give  the  distribution  function  N(w)  in  terms  of  the 
measured  current,  apart  from  a  multiplicative  constant  which  depends  on 
system  parameters,  for  example,  transparency  losses  at  grids,  detector 
sensitivity,  etc. 


We  also  wish  to  find  expressions  for  the  average  velocity  and  average  energy. 

In  getting  the  energy  we  must  remember  that  we  are  detecting  only  ions  whose 
y  and  z  components  of  velocity  are  essentially  zero.  We  assume  however,  that 
the  x  components  are  independent  of  the  y  and  z  components  as  regards  the 
random  part  of  the  distribution,  the  latter  being  superimposed  on  a  drift  velocity 
having  only  an  x  component.  Then  the  average  ion  energy  measured  may  be 
manipulated  as  follows: 


—  1 

E  *  —  m 


£  w2  N  (w)  d  w 
£  N  (w)  d  w 


E-  =  i.  m  Z[w  +  (w-w  )  3  2  N(w)  d1 


*  —  ~  TI*  2  +  2w  (w-w  )  +  (w-w)2]  N  (w)d1 


2  N 


1  -2,1  m 
-  m  w  +  -  — 

2  2  N 


£  (w  - 


—  2 

w  )  N  (w)  d  w 


where  N  =  ^N(w)dw 
^  N(w)  dw 


(19a) 


Thus  E  is  decomposed  into  a  center  of  mass  energy  and  the  mean  energy 
(x  component)  about  the  center  of  mass  motion.  According  to  the  assumption 
made,  the  y  and  z  components  of  motion  would  give  rise  to  energies  equal  to 
the  random  part  of  the  x  distribution.  Then  the  mean  total  energy  of  the  ions 
including  all  components  would  be 


®*total 


—  m  w  2  +  —  —  £  (w  -  w  )2  N  (w)  dw 

2  2  N 


(19b) 
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B.  FINITE  WIDTH  FOR  EMISSION  PULSE 


In  the  more  general  case,  where  there  is  not  a  one-to-one  correspondence 
between  escape  velocity  w  and  time  of  detection  t,  the  solution  (10)  cannot 
determine  the  distribution  of  ions  escaping  over  both  velocity  and  time  solely 
by  use  of  the  detected  current.  Additional  assumptions  concerning  the  ion 
distribution  must  be  made.  One  could  assume  certain  properties,  for 
example,  a  Maxwellian  distribution  of  initial  escape  velocities  corresponding 
to  some  temperature  which  may  or  may  not  vary  with  time  of  escape  or  one 
could  assume  an  escape  rate  which  could  be  a  square  pulse,  or  have  the  laser 
pulse  shape. 

Free  parameters  such  as  the  temperature  of  the  distribution  or  a  constant 
multiplicative  factor  of  the  number  escape  rate  could  then  be  adjusted  in 
accordance  with  the  detected  current  after  the  integration  called  for  in 
Eq.  (10)  has  been  performed. 

Perhaps  the  easiest  way  to  see  the  complexity  involved  in  using  the  solution 
(10)  when  the  ions  do  not  escape  the  plasma  in  an  essentially  zero  time 
interval  is  to  assume  that  they  escape  at  two  discrete  times  t'^  and  t 1 2 - 

2 

Let  fj  (t'j)  f  number  of  ions/cm  which  escape  at  t'j, 

2 

f2  (t!g)  =  number  of  ions /cm  which  escape  at  t'2, 

(w)  =  number/unit  escape  velocity  for  ions  escaping  at  t ‘ , 

N2  (w)  -  number/unit  escape  velocity  for  ions  escaping  at  V 


Consider  the  currents  measured  at  times  tj^  and  t2.  Let  w^,  \v9^, 

and  w22  be  the  escape  velocities  for  ions  emitted  at  t'^  and  t'2  and 
reaching  the  detector  at  t  and  t2<  The  second  subscript  refers  to  the  times 
t^  and  t.^,  and  the  first  subscript  to  the  limes  t'j  and  t'g.  Then, 


J(t^)  IdtJ  -  f^t'j)  N1  (wjj)  (dw^l  +  f2  (t'2)  N2  (w21)  |dw21| 

^2’  Kbh  (t*l»  N1  (w12Mdw12l  +  r2  (tV  N2  (w22)I  dw2ij 


(20) 


Now 


Now  consider  specific  cases  concerning  what  assumptions  need  be  made  to 
enable  Jvt;  to  determine  *he  remaining  quantities. 

(a)  Suppose  t  t’2,  N.(w),  N2(w),  are  known.  Compute  w  j  j,  \v;j, 

Wj2  and  w22  along  with  1 1  ,  etc,  from  (21)  and  (22).  Then  J(tj)  and 

1*7 

j<t2)  determine  fj  (t'j)  and  f2  (t'2>.  If  there  were  i  discrete  ion  pulses,  i 
measured  values  J  (t)  would  determine  f-  (t  .),  provided  the  N^(w)  were  all 
known.  More  experimental  values  would  over  determine  the  f.  (t'^)  and  a 
least  squares  fit  could  be  made. 

(b)  A  different  situation  holds  if  we  regard  t’j,  t'2,  f,  (t‘j)  ar.d  fg  (t'g)  as 
known.  Then  from  (20)  we  cannot  solve  for  Nj  (wjj),  Ng  (w2J),  Nj  (w12> 
and  N2(w22>.  Two  additional  relations  would  be  needed.  For  i  ion  pulses 
and  i  measured  values  J  (t.),  i  (i  -  1)  additional  relations  would  be  needed. 
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On  the  other  hand  if  we  assume  Nj(w),  say,  is  Maxwellian  (or  any  given 

distribution)  then  this  applies  to  N^w^),  N^Wjg)#  N^w^)  for  all  times 

t.  at  which  currents  J(t-)  detected  are  to  be  used.  If  we  assume  i  distri¬ 
bution  functions  N.(w)  corresponding  to  i  times  t'.  for  which  the  f.  (t’.) 

are  assumed  known  and  then  use  i  experimentally  determined  points  J(t-) 

we  would  be  able  to  determine  the  temperature  T,  in  each  of  the  i  distributions 
and  hence  in  the  limit  find  T(t’),  i.e. ,  the  temperature  of  the  assumed 
Maxwellian  distribution  as  a  function  of  time  of  escape  from  the  plasma. 

In  general,  one  would  not  trust  the  information  obtained  when  so  much  must 
be  assumed.  The  difficulty  is  that  one  detects  the  current  as  a  function  of 
only  one  variable,  time,  whereas  the  ions  escape  the  plasma  in  a  distribution 
of  two  variables,  velocity  and  time  (and  a  third,  space,  which  we  have 
ignored  because  the  ions  all  originate  or  escape  the  plasma  very  near  to  the 
target  surface;  a  small  potential  correction  could  be  applied  here  if  it  turns  out 
that  the  distance  from  the  target  is  at  all  significant). 

However,  while  more  must  be  assumed  to  exploit  the  general  solution  (10) 
than  must  be  assumed  for  a  delta  function  distribution  in  time,  the  delta  function 
is  a  more  drastic  assumption*  at  least  for  the  lighter  ions  where  the  time 
interval  of  the  mass  peak  is  not  greatly  longer  than  the  laser  pulse. 

Therefore,  the  most  reasonable  procedure  is  to  assume  a  delta  function  pulse 
in  time  as  a  first  approximation.  If  higher  accuracy  seems  achievable,  one 
could  then  use  the  general  solution  (10)  using  a  square  pulse  approximating 
the  laser  pulse,  or  a  triangular  pulse  to  approximate  it,  for  the  ion  emission 
per  unit  time  during  the  pulse  then  assume  a  Maxwellian  velocity  distribution 
with  a  temperature  varying  in  time..  T  h  e  r  s  p  p  e  a  r 

to  be  no  difficulty,  in  principle,  in  carrying  out  this  procedure  provided  a 
number  of  discrete  time  intervals  are  used  to  approximate  the  ion  emission. 

If  an  analytical  choice  of  pulse  shape  is  made,  ii  would  be  better  to  integrate 
(10)  at  the  outset  and  set  the  resulting  function  of  time  equal  to  J(t), 

However,  while  t  is  a  parameter  in  the  integral  in  equation  (10)  giving  the 
detected  current  as  a  function  of  time,  we  find  that  taking  the  simplest  case  ** 
a  square  emission  pulse  in  time  -- 

J(0,  w,  t)  *  o  (0,  w,  t)  w 
*  A{0f(w) 

with  Ait)  1  AQ  0<  t  <  T 

*  0  t  <  0 

*  >  *r 
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and  putting  this  into  the  solution  (10),  in  general,  one  cannot  perform  the 
integraMon.  Taking  f  fw)  to  be  a  (one  dimensional)  Maxwellian 


a  specia1  aMiough  perhaps  the  most  interesting  case,  does  not  help.  We 
arrive  at  an  expression 


T  a  function  cf  [  t  7  *w>  i  * 


Starting  wi’h  equation  «20)  Jft}  may  alternative1)'  be  expressed  as  an  Integra1 
ever  time  of  escape  from  the  plasma; 
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We  see  that  if  T  is  left  as  a  function  to  be  determined  by  J(t),  but  otherwise 
unspecified, 

T  =  T  [t  -  t(w)  ] 

or 

T  -  T  <t') 


the  integrations  cannot  be  carried  out.  If  T  is  taken  to  be  a  constant,  the 
integration  can  be  performed.  The  time  dependence  of  the  integral  is  due  to 
the  fact  that  the  limits  of  w  are  time  dependent,  not  to  anytime  dependence  of 
the  integrand.  (In  the  alternative  integral  over  t',  w  is  a  function  of  t  so  the 
integrand  hat  not  the  limits  are  time-dependent).  In  either  case  the  form  of 
the  function  J(t5  is  specified  by  the  integral  and  only  the  best  fit  of  the  detected 
current  J(t)  to  the  specified  form  can  be  obtained. 
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APPENDIX  II 


ION  EMISSION  FROM  LASER  IRRADIATED  TUNGSTEN 

(Paper  presented  at  the  F  irth  International  Quantum  Electronics  Conference 
Phoenix,  Arizona,  April  j66) 


In  this  paper  we  report  on  the  emission  of  positive  ions  from  tungsten  surfaces 
irradiated  with  the  output  of  a  Q-switched  laser  in  the  range  of  20-70  Mw/cm  . 
The  results  on  ion  emission  are  only  one  part  of  a  comprehensive  study  of  the 
effect  of  laser  irradiation  on  opaque  surfaces  at  power  levels  for  which  the 
surface  is  not  visibly  damaged.  Those  parts  of  the  study  dealing  with  cal¬ 
culations  of  surface  temperature  upon  irradiation,  as  well  as  electron  emission 
from  various  metallic  surfaces,  have  already  been  reported'-**  A  report  on 
the  neu’rai  particle  emission  is  now  being  prepared. 

Previous  workers  °  have  reported  the  observation  of  ions  of  the  irradiated 
metal  with  energies  up  to  a  1000  eVfor  light  fluxes  several  orders  of  magnitude 
higher  than  those  used  in  our  experiments.  At  those  power  levels  there  is 
considerable  evaporation  of  the  metal  and  most  of  the  emission  consists  of  ions 
of  the  substrate  material.  We  have  measured  the  energy,  masses  and  numbers 
of  ions  emitted  from  tungsten  under  conditions  where  most  of  the  emission  con¬ 
sists  of  species  adsorbed  on  the  surface.  We  find  most  of  the  ions  to  consist 
of  Na  and  some  of  which  have  energies  up  to  180  eV.  Other  ionic  species 
are  also  found.  The  ratio  of  the  ions  to  the  total  number  of  particles  indicates 
that  we  are  dealing  with  a  lightly  ionized  gas. 

Apparatus 

A  time -of -  flight  spectrometer  with  a  decelerating  grid  in  front  of  the  detector 
has  been  used  to  measure  the  charge-to-mass  ratio  as  well  as  the  energy  of  the 
ions  emitted  by  the  surface  during  irradiation.  The  instrument  is  very  well 
suited  tc  this  type  of  measurement  because  it  can  be  used  to  display  an  entire 
mass  spectrum  for  each  laser  pulse.  Such  a  feature  is  highly  desirable 
because  of  the  difficulty  in  obtaining  reproducible  quantitative  results  from 
successive  laser  pulses.  Figure  1  is  a  schematic  of  those  parts  of  the 
spectrometer  that  are  of  interest  here.  A  complete  description  of  the 
construction  details  has  been  given  elsewhere.  '•)  The  metal  to  be  studied, 
in  sheet  form,  is  used  as  an  electrode  and  held  1.  5  cm  away  from  a  one- 
me*er  long  drift  tube  terminated  in  a  disk  shaped  to  approximate  a  Rogowski 
surface.  On  the  opposite  end  of  the  drift  tube  is  a  decelerating  grid,  followed 
bv  a  DuMont  SPMO3-  301  electron  multiplier  detector. 
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INPUT 


Figure  2  Appendix  II  -  Experimental  Arrangement  for  Measurement  of  Time 

of  Flight 
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Figure  3a  Appendix  II  -  E  ■  e'  ■  +1000V.,  vertical  sensitivity  of  lower  trace  0.  2V/cm. 


Figure  3b  Appendix  II  -  E  ■  ■1000V.,  e'  ■  -USOV,  vertical  aenaitivtty  of 

lower  trace  0.05V/cm. 


In  operation,  the  target  surface  Is  kept  at  a  positive  potential  E  with 
respect  to  the  electrically  grounded  drift  tube.  A  positive  potential  E' 
is  also  applied  to  the  decelerating  grid,  while  the  front  of  the  detector  is 
at  a  negative  potential  E". 

Ions  produced  by  irradiation  of  the  target  with  the  output  of  the  laser  are 
accelerated  into  the  drift  tube,  where  mass  separation  is  achieved  simply 
because  ions  of  different  masses  have  acquired  different  velocities  under 
the  action  of  the  same  accelerating  potential  difference.  Upon  leaving  the 
drift  +ube,  the  ions  are  decelerated  by  the  positive  potential  applied  to  the 
grid,  and  if  they  surmount  the  potential  E1,  are  finally  accelerated  into  the 
detector.  The  multiplier  current  is  measured  across  a  220  ohm  resistor, 
giving  a  detector  time  constant  of  <  10  seconds. 

Ex  per  i  men t  al  Details 

A  ruby  laser  with  a  cr vptccyanine  Q-switch  was  used  in  all  the  experiments. 

The  output  of  the  laser  consists  of  a  nearly  triangular  pulse  with  a  full  width 
of  30  nanoseconds  at  the  base,  and  peak  power  of  6  Mw.  The  light  flux  at  the 
target  was  20  70  Mw/cm^  depending  on  the  lens  used  to  focus  the  beam. 

The  experimental  arrangement  is  shown  in  Figure  2.  The  laser  is  aimed  aT 
the  surface  through  a  sapphire  window  in  the  spectrometer  envelope,  at  an  angle 
of  53  degrees  ~o  the  normal.  The  signal  from  the  electron  multiplier  measuring 
the  ion  current  in  the  time  of --flight  spectrometer  is  fed  to  the  two  vertical 
inputs  of  a  Tektronix  555  dual  beam  oscilloscope.  Time  base  A  of  the  scope  is 
triggered  by  the  signal  from  a  Raytheon  917  phototube  monitoring  the  laser  pulse. 
T’tme  base  B  is  tr:ggered  internally  by  the  electron  multiplier  signal  after  a 
variable  delay  time  has  elapsed,  provided  that  B  has  been  preconditioned  by 
the  triggering  of  time  base  A.  The  gate  outputs  of  time  bases  A  and  B  are 
fed  into  the  start  and  stop  inputs  of  a  CMC-800A  digital  counter.  Different  mass 
pulses  can  be  measured  by  simply  varying  the  trigger  delay  on  time  base  B. 

The  output  of  the  917  photodiode  is  also  used  to  trigger  a  Tektronix  585 
oscilloscope  that  displays  the  output  of  a  Philco  4501  diode  monitoring  the 
laser  pulse. 

Ion  Energies 

Figure  3  shows  spectra  obtained  from  a  tungsten  surface  irradiated  with  a  flux 
of  50  Mw'cm  ,  and  held  at  a  positive  potential  E  =  +1000V.  The  upper  trace 
shews  the  laser  pulse:  the  lower  trace  is  the  multiplier  signal.  Time 
increases  to  the  left. 

Figuie  3a  was  taken  with  E'  =  E  -  -‘-1000V,  i.  e.  ,  no  net  decelerating  potential. 
Under  these  conditions  all  ions  emitted  by  the  surface  within  the  spectrometer 
aperture  ate  collected.  Figure  3b  is  taken  under  the  same  set  of  conditions 
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as  3a*  except  that  the  potential  of  the  decelerating  grid  has  been  increased 
to  E;  -  +1180V. ,  while  E  =  +1000V.  If  E'  is  further  increased  to  +1200V, 
while  holding  E  constant  at  +1000V,  all  peaks  disappear.  The  presence  of 
the  two  major  peaks  in  Figure  3b,  that  have  been  identified  as  Na+  and  K+, 
indicates  that  some  of  the  ions  have  energies  in  excess  of  180  eV.  On  the 
other  hand,  the  ability  to  block  the  pulses  from  reaching  the  detector  by 
increasing  E'  assures  that  they  are  charged  particles. 

There  are  two  other  factors  that  support  the  existence  of  high  energy  ions. 
First,  the  number  of  ions  per  peak  in  Figure  3b  is  only  1%  of  the 
corresponding  number  in  3a.  Second,  it  can  be  seen  that  the  narrowing  of  the 
peaks  of  Figure  3b  has  occurred  by  the  elimination  of  the  low  energy  tail 
in  the  peaks  of  3a. 

At  this  point  it  seems  worthwhile  to  consider  poss'ole  sources  of  spurious 
signals  that  might  appear  to  be  high  energy  ions  ir  the  spectrometer.  One 
may  rule  out  the  possibility  of  electrons  generated  anywhere  in  the  system 
arriving  at  the  detector,  since  they  would  have  to  surmount  a  potential 
barrier  of  3  5  kilovolts  immediately  in  front  of  the  multiplier.  Two  other 
possible  sources  of  spurious  signals  are  UV  emitted  by  the  heated  target,  and 
ions  generated  at  places  other  than  the  target.  The  former,  if  present,  would 
be  easily  identified  by  a  zero  time  of  flight  and  an  undiminished  signal  when  a 
strong  magnet  is  placed  near  the  target,  ions  generated  at  any  point  between 
the  target  and  decelerating  grid  will  not  have  the  full  energy  normally  given 
to  ions  in  the  accelerating  region.  Hence,  they  will  not  be  able  to  surmount 
the  decelerating  barrier,  since  E  >  E.  Ions  formed  past  the  decelerating 
electrode  will  be  accelerated  to  the  detector  and  recorded,  but  they  cannot 
be  stopped  by  the  decelerating  potential.  The  fact  that  we  are  able  tc  cut-off 
all  the  peaks  by  increasing  E  indicates  that  this  is  not  making  place. 

ion  Masses 


The  time  of- flight  0f  a  particle  of  mass  number  m,  from  target  to  detector, 
can  be  calculated  using  simple  equations  of  motion.  For  the  spectrometer 
dimensions  ^ the  ume  of  flight  for  a  singly  ionized  particle  is  given  by 


(1)  t  (usee)  -  1.44‘Vra 
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where  E  is  the  initial  energy  of  the  particles  (in  ev )  and  E,  Ef,  and  En 
are  the  absolute  values  of  the  voltages  'in  volts)  applied  to  the  electrodes, 
as  shown  in  Figure  2.  Figure  4  shows  the  time  of  flight  of  particles  of 
different  charge  to  mass  ratio  as  a  function  of  initial  energy  EQ,  for  E 
'  E  1000V  and  E1  -  -3000V.  It  can  be  seen  from  the  figure  that  the 
resolution  of  ’he  instrument  is  large  enough  that,  for  the  values  of  m/e  shown, 
no  serious  overlap  of  peaks  is  expected  for  initial  energies  in  excess  of  100  eV. 
Hence,  even  in  the  absence  of  an  energy  analyzer,  the  peaks  will  be  well 
enough  resolved  to  make  mass  assignments  possible.  This  fact  is  clearly 
demonstrated  by  the  well  resolved  peaks  shown  in  Figure  3., 

Table  1  compares  the  measured  and  calculated  times  of  flight  for  the  masses 
found  in  the  emission  from  the  tungsten  surface. 

TABLE  .1 

ION  MASSES  FROM  TIME  OF  FLIGHT 


Mass  ionic  Sp  ectes  Measured  TOF  (.us  Calculated  TOF  (us) 


•» 

H 

2  5 

2.  4 

7 

+ 

Li 

6.4 

6.  4 

12 

■4- 

C 

8.3 

8.  3 

18 

H2°+ 

10.  1 

10.2 

23 

Na  + 

11.  9 

11.6 

28 

co" 

13.  3 

12.  8 

39 

K 

15.  7 

15.0 

44 

co2‘ 

16.7 

16.0 

184 

wf 

33  0 

32.  7 

The 

sing 

fit*  a 

times  of  1 1  i gh *  were  cal 
:v  charged  thermal  ions 
.s  included' 

culated  using  *1 
.  A  corrcc’ion 

with  E  +10C0V,  and  assuming 
for  the  fringing  of  the  accelerating 

Tin  t-1  are  several  reasons  for  assuming  ’hat  the  ions  observed  are  singly 
char  god  Fast,  the  assumption  leads  one  to  mass  assignments  for  the 
obs <  rved  peaks  that  are  \n  agreement  with  the  resui’s  of  electron  and  laser  ^ 
In  tm  probing  of  polycrvs’a'iine  tungsten  surfaces  m  sector  and  quadrupole 
mass  spec*  i  cme*  ets  res  pec*  ive'y.  Second,  d  'he  partules  were  doubly  charged, 
one  wculd  expect  *o  see  a  aiger  peak  at  the  time  corresponding  to  the  arrival  of 
s mg.1  v  ehai  ged  ions  of  ’he  same  species  This  is  no»  found  in  any  of  our 
mi  is ut omen’s . 
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Figure  4  Appendix  II  *  Time  of  Flight  vs.  Initial  Ion  Energy 
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While  the  presence  of  large  quantities  of  alkali  metals  may  seem  surprising, 
they  are  a  common  contaminant  in  tungsten,  for  their  compounds  are  used  in 
the  processing  of  tungsten  ore  and  the  rolling  of  the  metal.  The  other  ionic 
species  found  are  known  to  adsorb  readily  or.  metal  surfaces. 

Number  of  Ions 

The  number  of  ions  of  different  species  that  are  emitted  per  laser  pulse  have 
been  measured  by  integration  over  peaks  such  as  those  in  Figure  3.  The 
calibration  procedure  for  the  spectrometer  hqs  been  given  elsewhere  + 
Typical  numbers  for  a  light  flux  of  70  Mw/cm2  are  9  x  10°  K*,  2  x  10®  Na 
and  2  x  10®  W  ions  respectively.  These  three  ionic  species  make  up  most 
of  the  emission  from  the  heated  surface.  The  relative  numbers  of  ions  do 
not  represent  the  relative  abundance  of  different  species  in  the  emission, 
since  the  alkali  metals  will  boil  off  the  tungsten  primarily  as  ions.  In  fact, 
measurements  of  the  neutral  particle  emission  made  in  a  quadrupole 
spectrometer'9)  for  the  same  laser  powers  indicate  that  there  are  approxi¬ 
mately  101®  particles  emitted  per  pulse.  This  also  means  that  we  are  dealing 
with  a  lightly  ionized  gas,  which  explains  why  the  ions  react  readily  to  external 
electric  fields. 


The  number  of  W  ions  measured  per  pulse  has  been  compared  with  the  number 
predicted  by  the  Saha  equation  ■  '  for  a  tungsten  surface  at  the  temperature  of 

the  target.  The  temperature  of  the  surface,  as  well  as  the  number  of 
electrons  ^mitted  at  that  temperature,  were  calculated  by  methods  previously 
described  .  The  number  of  tungsten  ions  was  then  calculated  using  known 
mass  evaporation  r  tes'11).  The  number  of  ions  predicted  by  the  Saha 
equation  using  the  values  obtained  from  the  above  calculations  is  10®.  Hence, 
we  see  that  the  numbers  of  ions  is  considerably  higher  than  would  be  expected 
from  a  purely  thermal  process. 


Inverse  Bremsstrahlung  Heating 

The  question  of  the  origin  of  the  high  energy  ions  and  their  acceleration 
mechanism  is  of  interest.  We  have  considered  a  number  of  possible  mechanisms 
To  date  there  is  no  satisfactory  description  of  the  processes  in  which  these  ions 
originate.  The  mechanism  which  produces  high  energy  ions  must  operate 
quite  generally,  since  they  are  produced  under  a  wide  variety  of  experimental 
conditions.  The  ion  velocities  obtained  in  the  present  experiment  are  in  fac* 
comparable^  to  those  observed  when  much  higher  laser  power  densities  are 
employed.  *  * 
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(12) 

Many  workers  have  considered  the  importance  of  absorption  by  inverse 
Bremsstrahlung  in  heating  a  laser  produced  plasma.  While  this  process  is 
capable  of  producing  considerable  heating  in  a  dense  plasma,  it  does  not  appear 
to  be  operative  in  producing  the  ions  observed  here^S),  However,  because  of 
the  great  present  interest  in  this  process,  both  for  producing  hot  dense  plasmas 
and  for  contributing  to  the  growth  of  laser  produced  gas  breakdown,  it  is  worth¬ 
while  to  develop  the  reasoning  leading  to  the  present  negative  result. 

The  early  portion  of  the  laser  pulse  vaporizes  material  which  is  emitted  as  a 
plasma.  Since  the  power  density  under  our  experimental  regime  is  relatively 
low  and  since  the  surface  is  not  extensively  damaged,  we  assume  the  material 
is  relatively  cool,  in  equilibrium  with  the  relatively  slowly  vaporizing  surface, 
and  is  only  slightly  ionized.  A  fraction  of  the  later  part  of  the  laser  pulse  is 
absorbed  in  the  plasma  by  absorption  of  a  photon  by  a  free  electron  in  the  plasma 
in  the  inverse  Bremsstrahlung  process.  The  previous  workers considering 
a  hot  dense  plasma  as  the  starting  point,  have  showed  that  this  interaction  can 
provide  considerable  heating,  to  the  point  where  hundreds  of  electron  volts  per 
ion  of  thermal  energy  are  produced. 

However,  these  results  are  not  necessarily  applicable  in  our  condi-ions.  We 
can  indeed  envision  a  heating  process  of  this  type  in  our  postulated  plasma. 

The  presence  of  a  small  amount  of  ionization  leads  to  a  small  amount  of 
heating  which  in  turn  leads  to  an  increase  in  ion  density  through  thermal 
ionization.  This  causes  an  increase  in  the  absorption  coefficient,  and  hence 
a  higher  rate  of  heating,  so  that  the  process  leads  to  a  rapid  runaway  heating. 
The  process  is  opposed  by  expansion  of  the  plasma,  which  lowers  the 
particle  densities,  and  by  radiation  of  energy  from  the  plasma,  which  lowers 
the  temperature.  The  ion  density  must  eventually  be  limited  by  the  fotal  number 
of  heavy  particles  present. 

The  absorption  coefficient  kv  for  absorption  of  light  of  frequency  v  by  inverse 
Bremsstrahlung  in  a  plasma  of  singly  changed  ions  of  density  nt  and  electrons  of 
density  ne  including  cooperative  effects  In  the  plasma  is  given  byU4) 

8rrn  n  e  l  n  A 

k  =  - - 

3  ev2  <2  nm  kT)3/  2 


where  e  and  m  are  the  electronic  mass  and  charge,  T  the  plasma  t^mpera’ut e, 
and  k  is  Boltzmann  s  constant.  Here  A  is  given  try 

jo.]3'2 

tt  n  / 

e  / 


3tt  n 

A  .  - 1 
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Over  a  broad  range  of  temperatures  and  electron  densities  covering  tne  range 
of  interest.  In  H  is  approximately  equal  to  10. 

The  expansion  of  the  plasma  is  assumed  to  proceed  as  a  complete  rare¬ 
faction  wave  We  assume  initial  conditions  for  the  material  with  respect 

to  particle  densities  and  ’emperature  We  assume  ’he  material  is  in 
equilibrium  wi*h  the  surface  and  fhat  the  ion  density  is  given  by  the  Saha 
equation’31^,  i.e.,  that  collisions]  processes  in  the  blcwcff  material  occur 
sufficiently  rapidly  to  keep  the  ion  concentration  in  equilibrium  with  the 
temperature  and  the  other  particle  densities-  This  assumption  is  open  to 
serious  question  for  our  experimental  situations  however,  if  it  does  not 
accuiateiy  describe  the  situation,  *he  actual  temperature  will  be  lower  than 
what  we  calculate  here  and  the  conclusions  of  this  analysis  will  be  unchanged. 

For  a  volume  of  blowoff  material  beginning  to  expand  near  the  target,  k  can  be 
calculated  and  thus  the  rate  of  heaung  due  to  the  absorption  of  radiation  can  be 
determined  from  the  initial  particle  densities.  Then  the  changes  in  particle 
densities  are  oalcula’ed  from  +he  equanons  describing  the  rarefaction  process 
and  the  Saha  relation.  The  hea’ing  rate  can  fhen  be  redetermined  from  fhe 
revised  value  of  k.  Thus  the  ionization,  expansion,  and  heating  of  a  volume 
element  of  blowoff  material  can  be  followed  as  a  function  of  time  starting  from 
its  initial  condition  at  the  surface  of  the  target  This  computational  routine 
has  been  carried  out  on  a  computer. 

Results  for  two  sets  of  initial  conditions  appear  in  Figure  5a  and  5b,  The 
initial  elec" ron  density  is  an  independent  parameter.  The  initial  neu  ral 
molecule  density  is  different  m  +he  two  pans  of  The  figure. 

We  see  tha*  for  considerable  heating  to  occur,  the  electron  density  must 
be  relatively  high.  A’  low  electron  der.sries  there  is  prac’icaily  no  heating 
effect.  As  electron  density  is  increased,  we  reach  a  point  where  a  rapid 
nonlinear  heating  occurs-  The  increases  in  temperature  and  ion  density 
interact  on  each  other  an  increase  in  either  quantity  drives  the  o’her  upwards 
even  faster.  There  appears  to  be  a  critical  electron  density  above  which  the 
process  is  possible  and  below  which  only  slight  heating  occurs. #  In  this  parti¬ 
cular  case,  an  initial  electron  density  of  the  order  of  1039  cm  ^  is  necessary 
to  rtaul*  m  this  process  occurring  during  a  30  nanosecond  laser  pulse. 

Changes  m  the  other  inrial  parameters,  the  'emperature  and  the  neu*ral 
par  Mole  density,  have  only  a  small  elf‘*c’  on  ’he  behavior  of  the  svstem, 
although  tlu-y  do  dc:»  imuu*  th**  value  of  tin  cr  deal  »  l*-<  tron  density.  Com* 
parson  of  Figure  5a  woh  Figure  5b  shows  *hut  the  critical  elec*  ran  density 
it.ct‘*ase*  as  initial  parttele  density  decreases  al*hough  the  varia*icn  is  not 
e<ttem«ly  sensitive  The  vVues  of  initial  mu»ra?  density  used  m  Figures  5a 
and 5b bracket  a  range  n  which  the  anus'  va’ue  mav  be  expect’d  to  fa  i  We 
shall  es'ima’e  an  upper  l*.mt»  to  this  va  =  ue  n**\* 
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Figure  5b  Appendix  U  -  Initial  Neutral  Particle  Denary  ’O19  cm 
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RATIO  OF  SON  OENSITY  WITH  ABSORPTION  TO  ION  DENSITY  WITHOUT  ABSORPTION 


We  may  calculate  an  upper  limit  to  the  initial  density  of  neutral  particles  ^ 
by  calculating  the  amount  of  material  vaporized  by  the  laser  theoretically 
Using  ‘echmaues  described  in  reference  i,  we  calculate  a  vaporized  depth  of 
5  x  10  7  cm  for  cur  presen*  experimental  situation.  At  a  thermal  velocity  of 
10b  cm  sec,  this  would  lead  to  a  neu’ral  particle  density’ of  4  x  10-  8  cm 
The  dep*h  of  material  removed  per  pulse  can  also  be  measured  by  striking  the 
same  spot  with  a  number  cf  ideirical  laser  pulses  and  then  dividing  the  depth 
of  the  hole  produced  by  the  number  of  pulses,  We  carried  ou*  ’his  measure 
ment  using  a  me’allograph.  Because  of  ack  of  sensitivity  in  measuring  very 
small  depths,  we  determined  only  an  upper  limrt  to  the*  depth  removed.  For 
the  presen*  experiment,  we  determine  tha*  the  dep*.h  removed  per  pu'se 
must  have  been  less  than  10  6  Cm,  ]n  reasonable  agreement  with  the  value 
calculated  above.  We  note  that  the  resu^s  in  Figures  5a  and  b  bracket  a 
range  of  ini’ial  neu*ral  particle  density  consis*ent  with  these  results. 

We  have  m  this  analys’s  shewn  *ha»  inverse  Bremss*  rahlung  can  lead  to 
hea*ing  of  the  blowoff  ma*eriai  under  certain  condnons,  the  most  important 
of  which  is  a  high  initial  elec*ron  density.  Hew*  ver.  for  our  e  <perimetr  all 
conditions,  the  required  electron  density  is  ’oo  irgh  *o  allow  for  'his  process 
to  proceed  From  measuremeirs  of ’he  thermionic  ems  sum  curren’  produceo 
by  a  Q  switched  laser wh  * 
the  surface,  we  esnma'e  a  i  ea 
cm  ^  for  our  experimental  condnons  The  numbers  >  f  electrons  ernned  wo-*' 
obtained  under  c end. ‘ions  m  winch  The  acce>ra* ing  voltage  was  large  »-nough  *o 
overcome  space  charge  effects. 

The  resul’.s  quoted  earlier  in  ‘he  section  on  numbers  of  icns  es*ub:’>h  ’ ha •  Hie 
gas  is  verv  liglvU  u  niz^d.  This  resu’y  taken  m  conjunction  wi’h  ’in  above 
upper  ';m!’  for  neutral  par’icle  emission,  supper*  *he  result  *ha*  the  initial 
elec*ron  density  mus*  be  orders  of  magnitude  lower  ’ban  lO^/cm^. 

Thus  the  pit'ure  most  consistent  wi’h  »  xpei  '.men*  a'  results  for  the  range  of 
laser  pewer  densities  20  70  megawa*ts  'cm 2  is  ’hat  of  a  slightly  ionized  gas 
The  ions  and  electrons  are  produced  along  wrh  a  larger  number  of  neutral 
a*oms  and  molecules,  The  in. nut  electron  densnes  are  too  low  ’o  lead  to 
plasma  lea*  mg  by  ’he  processes  described  above  The  caiculj*ed  rosuKs  do 
appear  ‘o  b-  qual oa’ i ve * \  s’milar  to  wha*  ’s  observed  with  higher  laser 
powers,  win  ><  ’lie  er  rival  particle  densities  may  be  attained.  Further 
Hives*-  gat  torts  of  *  lie  production  .me  chan’ sm  »n  •  his  rangt  o*  ia.-'u  powers  . 
should  apparently  consider  menus  of  producing  high  direct  d  uim  g*^  rutliei 
than  high  thermal  energies  lor  ’ire  >on> 
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FINAL  REPORT  CONTRACT  NO.  DA- 1 1 -022 -AMC- 1749(A)  MOD  NO.  2 
MECHANISMS  OF  LASER  SURFACE  INTERACTIONS 

Page  5:  10  microseconds /cm  should  read  8.  5  rnicroseconds/cm. 

Page  10:  Lines  10  through  20  should  be  replaced  by; 


where  R  is  the  radius  of  curvature  of  the  trajectory  in  cm,  E  is  the  particle 
energy  in  ev,  B  is  the  field  in  gauss,  and  m  is  the  mass  in  atomic  mass  units. 
For  B  -  500  gauss,  E  =  200  ev,  and  m  =  50  amu,  R  is  about  30  cm.  Since 
this  500  gauss  field  could  be  applied  over  an  area  encompassing  at  least  5  cm 
of  particle  path  length;  it  should  be  possible  to  deflect  such  ions  approximately 
0.  5  centimeter.  Since  this  distance  is  about  equal  to  the  separation  between 
the  quadropole  rods,  it  should  be  possible  to  deflect  these  particles  out  of 
the  system.  The  conclusion  is  even  stronger  for  particles  of  smaller  mass. 
Thus  it  should  be  possible  to  stop  even  relatively  fast  ions  with  the  magnetic 
field.  In  the  experiments  described  below  both  methods  (mass  filter  and 
magnetic  field)  were  used. 

Page  28:  Equation  11-11,  w  should  read  w* .  Third  line  from  bottom  should 
;  ou.d:  whore  vj'  -  initial  ion  velocity  plus  velocity  imparted  by  field  between 
target  and  grid. 

Pvigt  36:  w  should  be  replaced  by  w  in  both  the  ordinate  and  abscissa. 

Page  37:  w  should  be  replaced  by  w  in  both  the  ordinate  and  abscissa. 

6  5 

Page  42:  *0  x  10  in  the  abscissa  should  read  40  x  10  . 
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